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ABSTRACT 

Turbulent  flow  drag  reduction  in  recent  years 
has  received  considerable  attention  in  the  engineering 
literature.  While  many  authors  purport  to  have  quanti¬ 
tatively  defined  the  drag  reduction,  it  is  clear  that 
under  some  conditions  extremely  large  entry  lengths 
occur  and  therefore,  there  may  be  serious  errors  in 
their  measurements . 

Entry  lengths  and  contraction  losses  have  been 

determined  for  both  Newtonian  and  drag  reducing  non- 

Newtonian  fluids.  The  drag  reducing  fluids  considered 

were  0.01%  and  0.20%  by  weight  of  Separan  AP-30  in  water. 

Data  were  obtained  in  2  inch  and  1/2  inch  pipes  over 

3  5 

Reynolds  numbers  ranging  from  approximately  10  to  4  x  10  . 
Data  obtained  with  the  Newtonian  fluids  are  in  good  agree¬ 
ment  with  the  available  data. 

For  the  viscoelastic  materials,  the  contraction 
losses  owing  to  the  sharp-edged  6  to  1  cross-sectional 
area  reduction  were  found  to  be  as  much  as  75%  more  than 
Newtonian  fluids  at  similar  conditions.  The  measured  entry 
lengths,  although  suffering  from  an  appreciable  amount  of 
scatter,  have  been  found  to  be  much  larger  than  for  the 
turbulent  flow  of  Newtonian  fluids  and  pass  through  a 
maximum.  For  the  0.20%  solution,  entry  lengths  in  the  order 
of  200  diameters  were  found  at  a  Reynolds  number  of 
2  x  104 .  Similar  entry  lengths  were  found  for  the  0.01% 
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solution  at  a  Reynolds  number  of  3  x  10^.  These 
observations  have  important  consequences  when  it  is  noted 
that  the  experimental  pipelines  for  studying  these  mate¬ 
rials,  are  generally  designed  by  assuming  turbulent 
Newtonian  entry  lengths  of  about  50  diameters. 

Dimensional  considerations  indicate  that  the 
contraction  losses  depend  on  a  dimensionless  group  which 
includes  the  relaxation  time  of  the  material,  0.  It  was 
found  however,  that  this  single  dimensionless  group  was 
insufficient  to  interpret  the  results  qualitatively. 
Other  dimensionless  groups  such  as  the  generalized 
Reynolds  number  was  also  shown  to  be  insufficient. 
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CHAPTER  I 
INTRODUCTION 

1 . 1  GENERAL 

One  of  the  many  problems  confronting  the  engineer  is 
to  define  entry  length  and  contraction  losses  for  various 
flow  geometries  either  experimentally  or  theoretically. 
Recently,  owing  to  the  many  practical  implications,  a 
large  number  of  studies  have  considered  the  properties  of 
drag-reducing  additives  in  turbulent  flow.  In  this  turbulent 
regime  by  the  use  of  pressure  measurements  and  velocity  pro¬ 
files  determined  in  a  region-free  of  entry  effects  drag 
reduction  can  be  defined  experimentally. 

In  addition  to  drag  reduction,  analyses  as  discussed 
by  Philipoff  and  Gaskins  (60)  and  others  (28,49)  may  lead 
to  the  determination  of  the  normal  stress  difference.  The 
entry  region  and  associated  contraction  losses  are  closely 
coupled  with  the  study  of  heat  and  mass  transfer  (34)  in 
these  systems:  that  is,  heat  transfer  coefficients  have  been 
noted  to  be  lowered  in  proportion  to  the  amount  of  drag 
reduction. 

It  is  therefore  imperative  to  have  a  quantitative 
analysis  concerning  the  length  required  for  entry  effects 
to  become  negligible.  In  the  past  this  problem  has  been 
dealt  with  in  a  qualitative  sense  (48,  66)  indicating  entry 
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lengths  similar  to  those  in  the  Newtonian  transition  regime 
(63) .  No  quantitative  analyses  are  as  yet  available  to 
determine  experimentally  the  magnitude  of  the  entry  lengths 
relative  to  Newtonian  values. 

Therefore,  the  purpose  of  this  study  was  to  consider 
experimentally  the  turbulent  entry  lengths  for  drag- 
reducing  systems  in  round  tubes.  Equally  important  is  the 
magnitude  of  the  contraction  losses.  However,  before  such 
an  analysis  commences  an  understanding  of  the  analyses  of 
laminar  and  turbulent  flow  for  Newtonian  fluids  as  well 
as  the  analysis  for  laminar  non-Newtonian  flow  will  be 
discussed  in  some  detail.  This  background,  based  on 
available  analyses,  serves  as  a  framework  for  interpreting 
the  experimental  results  and  is  also  an  indication  of  the 
complexity  of  the  problem. 

1.2  ENTRANCE  LENGTHS  AND  CONTRACTION  LOSSES 

The  entry  length  may  be  defined  as  that  distance 
required  for  a  fluid  flowing  from  an  infinite  upstream 
reservoir  into  a  capillary  or  pipe  to  a  point  where  the 
downstream  flow  is  fully  developed.  Three  non-independent 
pieces  of  information  are  generally  sufficient  to  character¬ 
ize  the  flow  in  the  entry:  entry  lengths,  contraction  losses  and 
velocity  profiles.  As  indicated  on  Figure  (1.1)  for  laminar 
flow,  the  fluid  enters  the  pipe  from  a  large  reservoir  and 
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FIGURE  (1.1) 

NEWTONIAN  LAMINAR  VELOCITY 
PROFILE  DEVELOPMENT 
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a  uniform  velocity  distribution  is  generally  assumed 
to  prevail  at  the  start  of  the  pipe.  Due  to  the  fric¬ 
tional  resistance  at  the  wall  a  boundary  layer  is  created 
at  the  wall  and  gradually  encroaches  on  the  uniform  stream 
as  the  flow  proceeds  downstream.  Since  at  any  downstream 
position  the  total  mass  flux  remains  constant  the  flow  in 
the  undisturbed,  unsheared,  central  core  accelerates  to 
compensate  for  the  fluid  retardation  near  the  wall.  This 
change  of  the  steam  velocity  causes  a  greater  static  pres¬ 
sure  reduction  than  for  the  corresponding  well-developed 
pipe  flow.  Beyond  the  point  where  the  boundary  layer 
encompasses  the  cross  section  of  the  tube  the  velocity  profile 

continues  an  asymptotic  development. 

As  a  result,  entry  length  can  be  measured  as  a 
distance  downstream  where  the  velocity  profile  is  a  fixed 
percentage  of  being  fully-developed  (20,  33)  or  where  the 
pressure  gradient  becomes  a  constant  (30).  In  purely  vis¬ 
cous  non-Newtonian  fluids  the  velocity  profile  is  a  function 
of  the  flow  behavior  index  and  yields  as  a  result,  entry 
lengths  that  are  a  function  of  the  flow  behavior  index. 

Figure  (1.2)  indicates  qualitatively  the  entry  length 

(L  )  and  the  contraction  loss  (AP  )  for  a  hypothetical 
0  c 

run  when  pressure  drop  measurements  are  used  as  the 


criteria . 
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FIGUREQ.2) 

NEWTONIAN  PRESSURE  PROFILE 
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In  laminar  flow  of  Newtonian  fluids  the  entry 
lenth  may  be  expressed  as  a  function  of  the  Reynolds 
number  as  follows: 


L 

e 


D 


(1.2.1) 


in  which  K  is  a  constant  equal  to  0.0555  (20).  For 
purely  viscous  non-Newtonian  fluids  the  value  of  K  has 
been  reported  as  being  a  function  of  the  flow  behavior 
index  (6,  13,  21,  76).  For  viscoelastic  fluids,  entry 
lengths  have  been  observed  to  be  shorter  (28,  30),  or 
longer  (55,77)  than  those  for  Newtonian  fluids  at 
similar  Reynolds  numbers.  Metzner  and  White  (47) show  that 
for  viscoelastic  solutions  equation  (1.2.1)  can  be  re¬ 
written  for  flow  between  parallel  flat  plates  as: 


L  L 

(  ^  )  =  (  5^  )  +  a  (n,s)N  (1.2.2) 

b  VE  b  PV  ws 

where  the  second  term  represents  the  incremental  change 
in  entry  length  due  to  elastic  effects.  Depending  upon 
the  values  of  n  and  s  the  coefficient  a(n,s)  may  be  nega¬ 
tive,  zero,  or  positive.  The  variation  of  a(n,s)  thus 
offers  a  partial  explanation  why  entry  lengths  for  visco¬ 
elastic  solutions  in  laminar  flow  can  be  larger,  greater 
and  possibly  equal  to  their  counterpart  Newtonian  values. 
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Considering  Newtonian  laminar  flow  in  a  pipe 
free  of  end  effects  one  can  write  Poiseuille's  equation 
as  : 


AP 


DEV 


pv2/g. 


=  32.0 


RN 


Re 


(1.2.3) 


The  above  can  be  extended  to  purely  viscous 
non-Newtonian  fluids  by  replacing  the  Reynolds  number 
by  the  generalized  Reynolds  number  (46)  as  indicated 
by  Drexler  (28).  The  excess  pressure  arising  from 
several  effects  can  be  incorporated  into  equation 
(1.2.3)  as  follows  (28): 


AP 


TOT 


pV2/g, 


=  32.0 


+  <j> 


RN 


Re 


(1.2.4) 


Here  APTQT  is  the  total  pressure  drop  over  the  length 
L  and  <j>  is  referred  to  as  the  Hagenbach  correction 
factor.  With  reference  to  equation  (1.2.4)  the  contrac¬ 
tion  loss  can  be  empirically  expressed  as: 


(1.2.5) 


■ 
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Astarita  and  Greco  ( 6  )  recently  have  expressed 
a  similar  relationship  between  the  dimensionless  contrac¬ 
tion  loss  and  the  Reynolds  number  as: 

AP  Kn 

c  1 

- a -  =  -  +  K„  (1.2.6) 

pv  /gc 

K1  . 

where  - —  is  the  Couette  correction  factor  and  K9  is 
equivalent  to  <p ,  the  Hagenbach  correction  factor.  Their 
data  and  analysis  indicate  that  for  any  Reynolds  number 
less  than  a  critical  Reynolds  number  of  about  146  the 
dimensionless  contraction  loss  becomes  inversely  propor¬ 
tional  to  the  Reynolds  number.  The  proportionality  constant 
was  experimentally  determined  to  be  795.  Beyond  the 
critical  Reynolds  number  the  dimensionless  contraction 
loss  was  found  to  be  a  constant  equal  to  5.48.  Both 
and  were  reported  for  a  specific  ratio  of  reser¬ 
voir  to  tube  diameter  3,  and  one  would  assume  that  these 
variables  could  be  highly  sensitive  to  the  geometry  of 
the  contraction  section.  Other  authors  have  reported 
the  Hagenbach  factor  as  a  constant  as  indicated  by 
Table  (1.1).  It  is  interesting  to  note  that  these  semi- 
empirical  values  grossly  underestimate  the  value  of 

contraction  loss.  The  semi-empirical  approaches 
assumed  that  at  the  pipe  entrance  or  at  some 
imaginary  section  upstream  there  exists  a  flat 
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TABLE  (1.1) 

*  . 

Newtonian  Hagenbach  Correction  Factors 


Author 

Astarita 

Bogue 

Campbell  and  Slattery 
Christiansen  and  Lemmon 
Collins  and  Schowalter 
Tomita 


Reference 

6 

6 

5.48 

13 

2.16 

18 

2.18 

20 

2.274 

21 

2.33 

75 

2.20 

.  , 


- 
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velocity  profile.  In  reality,  the  experimental  evidence 
of  Uebler  (87)  and  the  analysis  of  Wang  and  Longwell  (78) 
show  that  the  velocity  profile  is  not  flat  for  a  sharp- 
edged  contraction.  In  addition,  the  presence  of  a  vena 
contracta  (6,  58)  downstream  of  the  contraction  would  prove 
to  be  impossible  if  the  profile  was  indeed  flat. 

The  spectrum  of  values  of  $  is  due  in  part  to 
the  other  assumptions  and  the  particular  analysis;  for 
example,  solving  the  equations  of  motion  by  neglecting 
no  terms  (78),  by  Pohlhau sen's  Momentum  Integral  Technique 
(13) ,  by  the  introduction  of  a  stream  function  (21)  or 
a  similarity  transformation  (1) .  Significantly,  the  usual 
assumptions  that  the  hydrostatic  pressure  is  a  function 
of  axial  position  only  and  that  radial  momentum  transfer 
by  convection  is  negligible,  have  been  the  cause  of  some 
controversies  (23,24,78). 

Extending  the  analysis  of  laminar  flow  contraction 
losses  to  non-Newtonian  fluids,  Astarita  (7  )  determines 
that  for  both  weakly  and  highly  elastic  solutions  the 
dimensionless  contraction  loss  was  inversely  proportional 
to  the  Reynolds  number  (as  indicated  by  equation  1.2.6). 

The  constant  of  proportionality  was  suggested  to  be  a 
function  of  the  flow  behavior  index  only,  and  have  values 
smaller  and  larger  than  Newtonian  fluids.  In  turbulent 
flow  there  are  few  analyses  available  which  provide  a 
general  prediction  of  entry  lengths  and  the  attendant 
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FIGUREdo  3) 

NEWTONIAN  TURBULENT  VELOCITY 
PROFILE  DEVELOPMENT 
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contraction  losses. 

Figure  (1.3)  indicates  qualitatively  the  features 

of  the  entry  in  the  turbulent  regime.  In  this  regime, 

the  boundary  layer  remains  laminar  only  up  to  a  transition 

point  which  occurs  at  a  boundary  layer  Reynolds  number 

5 

roughly  equal  to  10  .  Beyond  this,  the  boundary  layer 
is  turbulent.  A  large  portion  of  the  excess  pressure 
loss  occurs  near  the  pipe  entry  and  how  this  pressure 
loss  is  distributed  between  the  two  regimes  is  of  yet 
not  clearly  defined.  However,  with  accurate,  closely 
spaced  pressure  readings  near  the  entry,  one  can  in 
principle  predict  the  downstream  location  where  transi¬ 
tion  occurs.  Olson  and  Sparrow  (59)  use  this  device  and 
show  as  a  result  that  a  vena  contracta  is  present  for  a 
sharp-edged  contraction  while  for  a  smoothed  entrance  the 
laminar  portion  was  observed  to  extend  a  considerable 
distance  downstream.  Experimental  evidence  indicates  that 
entry  lengths  of  25  and  15  diameters  can  result  from  non- 
induced  and  induced  turbulence  in  the  entry  geometry  for 
a  N  =  5.0  x  105. 

JA0 

Latzko  (35)  determined  that  the  entry  length  for 
induced  turbulence  could  be  represented  as: 

Le 

—  =  0.693  N  1/4  (1.2.7) 

D  Re 


: 
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However,  the  above  formula  yields  shorter  entry  lengths 
that  have  been  found  experimentally.  Barben  and  Jones 
(11) determined  that  for  a  NRe  =  4.0  x  105  the  velocity 

profile  was  not  developed  at  41  diameters  downstream  while 
Latzko's  theoretical  formula  would  have  given  an  entry 
length  of  17  diameters.  Kirsten  (38)  measured  entry  lengths 
of  50  to  100  diameters  depending  upon  the  Reynolds  numbers 
while  Nikuradse  (35)  obtained  fully-developed  flow  at  25 
to  40  diameters  for  induced  turbulent  flow;  the  40  diameter 

5 

value  was  obtained  at  a  N  =  9.0x10. 

A  more  recent  analysis  was  made  by  Ross  (62)  who 
predicted  pressure  losses  for  Newtonian  fluids  by  consider¬ 
ing  the  integral  momentum  and  continuity  equations  and 
utilizing  available  boundary  layer  equations.  The  results 
indicated  that  the  pressure  losses  in  the  entry  were  a 
weak  function  of  the  Reynold  number  for  Reynolds  numbers 
in  the  order  of  10^.  Theoretically,  this  indicates  that 
the  dimensionless  contraction  losses  are  constant. 

Generally,  the  available  data  for  turbulent  flow 
indicate  that  the  entry  lengths  vary  widely  and  in  a 
complex  way  on  the  geometry  of  the  entrance,  the  Reynolds 
number  and  even  on  the  means  of  measurement.  An  additional 
although  dependent  factor,  which  has  been  largely  ignored 
in  the  papers  mentioned  previously,  is  the  turbulence  level 
of  the  entering  stream.  The  importance  of  turbulence  in 
the  free  stream  as  well  as  in  the  entrance  geometry  is 
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concisely  discussed  in  a  recent  paper  by  Dalla  Lana  and 
Christiansen  ( 22) . 

At  pipe  Reynolds  numbers  near  2100,  or  equivalently 
"transition  Reynolds  numbers",  entry  lengths  larger  than 
those  at  high  Reynolds  numbers  have  been  observed  (  63)  . 
These  data  are  also  presented  in  condensed  form  in  Hinze 
(35).  In  the  present  context,  these  data  may  be  especially 
important  as  some  evidences  exist  to  show  that  for  the  flow 
of  viscoelastic  fluids,  transitional  flow  may  prevail  to 
Reynolds  numbers  in  excess  of  30,000  (66). 
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CHAPTER  II 

BACKGROUND  AND  THEORY  FOR  DRAG-REDUCING  FLUIDS 
2 . 1  GENERAL 

Unlike  Newtonian  fluids  which  are  capable  of 
exhibiting  only  viscous  responses,  a  viscoelastic  solu¬ 
tion  is  capable  of  exhibiting  both  viscous  and  elastic¬ 
like  responses  to  a  deformation.  These  responses  depend 
in  general  upon  the  instantaneous  deformation  rate  and  the 
total  deformation.  It  has  been  shown  in  recent  literature 
(5,  50,  51)  that  the  elastic  nature  may  be  considered  to 
be  governed  by  what  has  been  termed  a  Deborah  and  a 
Weissenberg  number  (3,  35,  67).  The  Deborah  number 
representing  the  ratio  of  the  natural  time  of  the  fluid 
to  the  residence  time  of  the  material  in  the  process 
defines  whether  the  flow  is  one  in  which  the  elastic 
properties  of  the  fluid  are  likely  to  predominate.  For 
rapid  deformations:  that  is,  where  the  residence  time 
is  much  smaller  than  the  natural  time  of  the  fluid  the 
elastic  properties  of  the  solution  predominate.  At  the 
other  extreme,  the  fluid  relaxes  relatively  quickly  as 
it  moves  through  the  flow  field  so  as  to  adjust  substan¬ 
tially  to  the  immediate  surroundings.  In  this  limiting 
case,  the  fluid  becomes  indistinguishable  from  a  purely 
viscous  fluid.  The  Weissenberg  number  representing  the 
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intuition  as  a  parameter  in  addition  to  the  Reynolds 
number  defined  as  the  ratio  of  inertial  to  viscous  forces . 

In  the  immediate  vicinity  of  the  sudden  contraction/ 
and  within  the  developing  boundary  layer  the  rate  of  defor— 
mation  experienced  by  a  fluid  element/  changes  rapidly  with 
position  downstream.  Equivalently,  from  a  Lagrangian  point 
of  view /  the  rapid  changes  in  deformation  rate  may  be  viewed 
as  rapid  changes  with  respect  to  time.  Thus,  as  has  been 
shown  in  an  earlier  boundary  layer  analysis  (44) ,  large 
values  of  the  Deborah  number  may  be  expected  to  occur  at 
the  leading  edge  of  the  tube  and  for  some  distance  downstream 
of  the  leading  edge.  Furthermore,  it  has  been  shown  that 
the  large  Deborah  numbers  may  lead  to  a  substantial  thicken¬ 
ing  of  the  boundary  layer  (44,75,80)  which  in  the  case  of  flow 
in  the  entrance  of  a  tube  or  channel  (47,51)  may  be  expected  to 
result  in  a  solid-like  region  immediately  inside  the 
entrance . 

Similar  effects  owing  to  elasticity  occur  in  the 
central  core  where  continuity  requires  that  acceleration 
must  occur  as  governed  by  the  growing  boundary  layer. 

The  deformations  in  this  region  have  been  termed  elonga- 
tional  flow  (5,77)  where  it  has  been  shown  as  in  the  previous 
paragraph  that  the  limiting  or  high  Deborah  number  behavior* 
of  the  fluid  is  a  solid-like  response.  For  the  viscoelastic 
fluids,  for  which  a  thickened  boundary  layer  is  expected 
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to  occur,  there  will  be  a  tendency  to  increase  acceleration 
of  the  central  core.  At  sufficiently  high  Deborah  numbers 
within  the  core,  the  acceleration  would  be  inhibited  and  the 
growth  of  the  boundary  layer  must  be  modified  in  order  to 
accommodate  the  resistance  to  high  deformation  rates. 

The  coupled  and  competing  effects  of  a  tendency 
toward  a  thickened  boundary  layer  and  a  resistance  to  large 
accelerations  within  the  core  region  imply  that  large  changes 
may  occur  in  the  pressure  profile  of  viscoelastic  fluids  as 
compared  to  those  of  Newtonian  fluids.  The  presence  of  a 
solid-like  region  implies  increased  contraction  losses  must 
occur  for  viscoelastic  fluids.  Both  effects  are  dependent 
on  a  Deborah  number  of  the  flow  which  in  turn  depends 
directly  on  the  fluid  relaxation  time. 

2.2  CONSTITUTIVE  EQUATIONS 

Before  analyzing  entry  length,  a  constitutive 
equation  must  be  chosen  to  represent  the  flow  behavior. 

In  the  available  literature  Metzner  and  Reed  (26,46) 
develop  and  employ  a  power  law  model  for  purely  viscous 
solutions.  This  power  law  is  expressed  as: 

8V  n' 

T  =  K'  (  —  )  (2.2.1) 

W  D 

where  n'  is  the  flow  behavior  index  and  K'  is  the 
consistency  index.  Qualitatively,  the  larger  the  value 
of  K'  the  more  "viscous"  is  the  solution  while  n' 
describes  the  degree  of  non-Newtonian  behavior.  Using 
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equation  (2.2.1)  and  the  Fanning  friction  factor  defi¬ 
nition  the  generalized  Reynolds  number  arises  naturally 


Dn'v2  n'p 


g  K  '  8  n  '  ~ 1 


(2.2.2) 


For  this  study  the  fluid  properties  n'  and  K'  have 
been  chosen  to  adequately  represent  the  viscous  nature  of 
the  polymeric  solutions.  Many  constitutive  equations  have 
been  developed  to  portray  the  elastic  properties  of  inter¬ 
est  and  two  of  these  constitutive  equations  discussed  in 
the  literature  are  the  Coleman  and  Noll  and  the  Convected 
Maxwell  models. 

The  Coleman  and  Noll  constitutive  equations  as 
discussed  in  the  literature  (81,82)  have  stresses  not 
dependent  on  the  instantaneous  deformations  but  on  the 
entire  history:  that  is,  strains  in  the  distant  past  have 
lesser  effects  than  in  the  recent  past.  The  mathematical 
elegance  of  this  formulation  however,  is  immediately  lost 
when  applying  it  to  problems  of  interest.  The  truncated 
constitutive  equation  such  as  the  second  order  approxi¬ 
mation  is  only  capable  of  predicting  normal  stress  behav¬ 
ior  at  deformation  rates  that  are  so  low  as  to  be  of 
little  importance  in  engineering  problems.  Furthermore, 
a  constant  viscosity  coefficient  is  maintained  thus,  limiting 
its  usefulness  to  only  slightly  non-Newtonian  fluids.  In 


•  e.oi  t:  -j  .  :  o- .  -  ...  -2  >  ti99  [ 


' 


19 


comparison,  the  third  order  approximation  introduces  a 
variable  viscosity  but  does  not  solve  the  problem  insofar 
as  normal  stresses  are  concerned.  Higher  order  approxima¬ 
tions  of  fourth  order  and  above  only  add  to  the  complexity 
of  the  problem  due  to  the  presence  of  a  larger  number  of 
material  coefficients  that  must  be  determined  experimentally. 
For  the  problem  at  hand  the  limited  range  of  applicability 
of  the  Coleman  and  Noll  constitutive  equation  suggests 
another  equation  must  be  chosen  to  portray  the  fluid 
properties . 

Metzner,  White  and  Denn  (51)  discuss  in  detail 
the  Maxwell  model  and  conclude  that  it  is  capable  of 
predicting  all  experimental  effects  at  least  qualitatively, 
except  for  the  assumption  of  a  zero  second  normal  stress 
difference.  Recently,  various  authors  (37,71,72)  show 
that  in  reality  this  second  normal  stress  difference  is 
not  zero  and  seems  to  be  related  to  the  first  normal  differ¬ 
ence  by  a  constant  that  is  yet  not  clearly  defined.  Measure¬ 
ment  of  the  first  normal  stress  difference  generally  employs 
a  jet  thrust  technique  (48,83)  for  high  shear  rates  and  a  cone 
and  plate  viscometer  (86)  at  low  shear  rates.  Middleman 
and  Gavis  (53)  point  out  drawbacks  to  the  former  procedure. 
Other  constitutive  equations  of  equal  complexity  have  been 
presented  however,  no  single  constitutive  equation  of  a 
tractable  complexity  is  available  to  universally  depict 
the  properties  of  polymer  solutions. 
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As  a  first  approximation  it  is  assumed  that  the 
convected  Maxwell  model  can  portray  the  fluid  properties 
such  that 


IP.  . 

9  C  13  =  2pdi  (2.2.3) 

P5  .  .  +  T.  . 

13  13 

Here,  P  is  an  arbitrary  hydrostatic  pressure  and  - —  is 
the  contravanant  convective  derivative  (51).  Thus, 
as  a  first  approximation  the  significant  elastic  proper¬ 
ties  can  be  assumed  to  be  governed  by  the  relaxation 
time  0 . 


where 


P.  . 
ID 


P.  . 
ID 


2.3  DRAG  REDUCTION  AND  VELOCITY  PROFILES 

Velocity  profiles  for  polymeric  solutions  in 
turbulent  flow  can  be  used  as  a  criteria  in  the  determina¬ 
tion  of  entry  lengths  and  in  determining  the  intermittency 
of  the  flow  field.  Recently,  in  developed  flow:  that  is, 
in  a  flow  field  assumed  to  be  free  of  all  entry  and  exit 
effects  experimental  data  indicate  that  the  velocity 
profile  can  be  expressed  as  (16,17,69): 


u+  =  a  log  y+  +  B  (  — —  )  -  C(C,f) 


v 


(2.3.1) 
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The  important  features  of  the  velocity  profile 
equation  are  discussed  in  detail  by  Ernst  (29)  ,  and 
other  authors  (52,69).  The  importance  of  this  equation 
is  that  coupled  with  the  definitions  of  the  Fanning 
friction  factor  and  average  velocity  a  unique  friction 
factor-Reynolds  number  correlation  can  be  developed 
(32,41,42,69).  This  relationship  quantitatively  por¬ 
trays  the  ability  of  the  polymeric  solutions  to  reduce 
frictional  losses  (19,31,66,70).  The  analysis  of  this 
phenomena  has  been  documented  experimentally  by  several 
authors  (3,45,67)  where  a  significant  reduction  of  the 
turbulent  drag  as  compared  to  the  flow  of  a  Newtonian 
fluid  under  similar  conditions,  is  noted.  Qualitatively, 
the  drag  reduction  is  displayed  as  a  friction  factor- 
Reynolds  number  plot  having  values  of  the  friction 
factor  as  much  as  an  order  of  magnitude  lower  than  corres¬ 
ponding  Newtonian  values  and  in  addition  exhibit  a  strong 
diameter  dependence. 

Recent  measurements  indicate  that  drag  reduction 
may  be  considered  to  be  a  function  of  elastic  to  viscous 
forces  (45)  .  Some  investigators  suggest  that  a  stabili¬ 
zation  of  the  flow  to  Reynolds  numbers  in  excess  of  2100 
results  directly  from  the  fluids  elastic  properties.  Not 
all  polymers  are  drag-reducing  (65).  Some  polymers,  such 
as  the  100  ppm  considered  in  this  study  which  has  a  flow 
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behavior  index  of  unity  (58) ,  yield  significant  drag 
reduction.  Also,  data  based  on  concentrations  as  low 
as  10  ppm  (58)  indicate  that  significant  drag  reduction 
is  present  even  for  minute  quantities  of  polymer. 

2.4  PRESSURE  LOSSES  AND  DIMENSIONAL  ANALYSIS 

The  measured  presssure  drop  from  an  infinite 
upstream  reservoir  to  a  point  downstream  is  composed  of 
frictional  losses  in  the  upstream  reservoir,  in  the  entry, 
a  loss  due  to  the  developing  profile  and  frictional  losses 
in  the  fully-developed  section.  In  laminar  flow,  these 
effects  as  well  as  elastic  contributions  can  be  summarized 
as  ( 56 ) : 


AP 


TOT 


A (KE)  +  FUpStream  +  Fdev  +  Fentry  + 


Fe+AE 


(2.4.1) 


In  the  above, F_  +  AE  are  contributions  owing  to 

E 

~the  elastic  properties  of  the  fluid.  The  term  E  is  the 
internal  elastic  stored  energy  per  unit  mass  and  F£  is 
the  dissipative  loss  due  to  any  unusual  flow  patterns 
caused  by  elasticity. 

In  viscoelastic  fluids  the  unique  ability  to 
store  and  release  quantities  of  energy  is  one  phenomenon 
yet  to  be  explained  in  detail  and  Astarita  (4)  suggests 
that  an  a  priori  understanding  of  the  thermodynamics  of 
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these  systems  is  necessary.  However,  it  is  not  clear 
that  these  elastic  effects  are  additive.  There  exists 
the  possibility  of  describing  the  total  pressure  loss  as 
a  function  of  the  non-elastic  terms  where  the  function 
might  depend  upon  the  fluid  relaxation  time  9.  The 
non-elastic  terms  in  equation  (2.4.1)  represent  the 
kinetic  energy  difference  for  a  fluid  that  flows  from  an 
infinite  reservoir  to  some  point  downstream  in  fully- 
developed  flow.  F  upstream  represents  the  frictional 
losses  in  the  upstream  reservoir  that  are  for  most  cases 
negligible  when  compared  to  Fdgv  and  Fentry:  that  is,  the 
frictional  losses  in  developed  flow  and  in  the  entry 
region. 

Feig  (  30)  and  Drexler  ( 28)  try  to  theoretically 
express  the  elastic  stored  energy,  E,  as: 

E  =  P:S  =  Z  2  P.,S..  (2.4.2) 

i  j  J  J 

where  P. .  is  defined  in  the  normal  manner  and  S. .  is  a 
13  13 

strain  tensor  that  results  from  a  deformation  in  a 
perfectly  elastic  solid  subjected  to  simple  shear  (  39)  . 

In  the  immediate  vicinity  of  the  contraction  this  assump¬ 
tion  may  be  valid  due  to  the  solid-like  behavior,  however , 
further  downstream  the  flow  field  is  free  from  any  solid¬ 
like  behavior  and  in  this  region  equation  (2.4.2)  is  no 
longer  applicable. 
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In  turbulent  flow,  while  in  principle  the  same 
considerations  apply  as  for  laminar  flow,  the  detailed 
analyses  of  entrance  phenomena  present  difficulties  that 
have  not  been  overcome  even  for  laminar  flow.  For 
viscoelastic  fluids  the  problem  is  of  course  complicated 
by  the  complex  nature  of  the  fluid,  thus,  dimensional 
analysis  will  be  applied  as  an  aid  to  interpreting  the 
experimental  results . 

Assuming  that  the  contraction  loss  is  a  function 
of  the  usual  variables  in  addition  to  the  flow  behavior 
and  consistency  indices,  and  the  relaxation  time  we  have 


AP  =  A(D,V,p,K' ,n' ,6,g  ) 
c  ^ 


(2.4.3) 


and  similarly, 


Lg  —  A1 (D ,V, p ,K 1 ,n ' , 0 ,gc) 


(2.4.4) 


Applying  the  theorem  we  obtain 


AP 


pv2/g. 


=  A 


g  K' 


ve 


DnV-n'p  ’  D  ' 


n 


(2.4.5) 


and 


=  A- 


D 


g  K’ 


Dn,V2  np  D 


ve 

vJ 

f  t  ^ 


(2.4.6) 


In  general,  these  dimensionless  groups  can  be 
rearranged  in  any  manner  provided  their  total  number  remains 
unchanged.  As  presented  the  three  dimensionless  groups 
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the  Reynolds  number,  a  type  of  Deborah  number  and  the 
flow  behavior  index. 

Since  Newtonian  fluids  possess  small  relaxation 
times  the  only  important  parameter  is  the  Reynolds  number. 

It  has  been  shown  that  the  dimensionless  contraction  loss 
can  be  either  inversely  porportional  to  or  independent  of 
the  Reynolds  number  (6).  The  dimensionless  entry  length 
presented  in  Chapter  I  indicates  a  direct  relationship 
with  the  Reynolds  number. 

For  viscoelastic  solutions  the  importance  of  any 
one  of  the  three  dimensionless  groups  is  not  known.  Through 
the  use  of  various  polymeric  solutions  possessing  a  range 
of  elastic  properties  one  can  determine  experimentally, 
relationships  between  the  dimensionless  contraction  loss 
or  entry  length  and  the  dimensionless  groups.  The  relaxa¬ 
tion  time  used  in  the  definition  of  the  Deborah  number  is 
that  which  is  determined  in  SLSF.1  Using  equation  (2.2.3) 
the  relaxation  time  in  this  case  is  defined  as: 


(P 


0  = 


11 


P22}w 


2T  )  (2.4.7) 

w  D 

Qualitatively,  the  relaxation  time  determines 
how  rapidly  the  stresses  will  decay  to  zero  after 
the  specimen  is  suddenly  stretched  and  held  at  a 
constant  deformation.  It  then  becomes  apparent 


1 


SLSF  refers  to  steady  laminar  shear  flow 
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that  the  flow  field  has  a  spectrum  of  values.  Near  the 
entry  as  previously  discussed  the  rapid  changes  in  defor¬ 
mation  may  be  viewed,  from  a  Lagrangian  viewpoint,  as  rapid 
changes  with  respect  to  time.  In  this  limit,  the  calcula¬ 
tion  of  relaxation  time  using  equation  (2.2.3)  becomes  complex 
since  it  now  becomes  a  function  of  time.  However,  these 
relaxation  times  are  related  to  those  defined  by  equation 
(2.4.7)  which  can  be  evaluated  experimentally.  The  Deborah 
number  will  thus  be  defined  using  the  relaxation  times 


obtained  in  SLSF. 
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CHAPTER  III 

3.1  EXPERIMENTAL  EQUIPMENT  AND  PROCEDURE 

3.1.1  GENERAL 

The  experimental  program  consisted  of  measuring 
pressure  losses  and  entry  lengths  for  dilute  aqueous 
solutions  of  a  polymer  additive  Separan  AP-30,  as  it 
flowed  through  a  sudden  6  to  1  diameter  contraction. 

Pressure  taps  were  located  before  and  after  the  con¬ 
traction.  Figure  (3.1)  is  a  schematic  representation 
of  the  experimental  equipment. 

3.1.2  PUMP  AND  MEASURING  SYSTEM 

The  pump  used  in  this  study  was  a  Moyno  2L10H 
having  a  variable  speed  drive  capable  of  delivering  a 
190  Imp. Gal ./Min.  output  flow.  To  prevent  the  discharge 
pressure  from  fluctuating,  a  surge  system  consisting 
of  two  10  foot  schedule  80  P.V.C.  pipes  were  connected 
to  the  discharge  line.  This  limited  the  output  pressure 
fluctuations  to  a  ±5%  deviation. 

Flow  rates  were  measured  using  a  2"  Foxboro 
Magnetic  Transmitter  in  a  vertical  position  and  a 
Model  696  Magnetic  to  Current  Converter.  These  were  recorded 
using  a  Beckman  Type  S-II  Dynograph.  For  recording 
purposes,  the  rated  10  to  50  milliamp  output  signal  was 
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converted  to  a  1  to  5  volt  signal  by  connecting  a  100 
ohm,  2%  resistor  across  the  output  leads.  Calibration 
consisted  of  employing  a  weigh  tank  for  small  flow  rates 
up  to  65  Imp . Gal . /Min .  and  friction  factor  measurements 
for  high  flow  rates  ranging  from  30  to  190  Imp .Gal ./Min. 
These  measurements  coupled  with  the  well-known  Von  Karman 
equation  determined  the  high  flow  rate  calibration.  Data 
for  the  flow  meter  calibration  are  presented  in  Appendix 
B,  Tables  1-3.  As  indicated  in  the  Appendix,  the  maximum 
deviation  in  flow  rates  is  within  1^%  of  the  maximum 
output . 

3.1.3  ENTRANCE  REGION  AND  PRESSURE  TAPS 

For  the  2"  test  section  the  entry  region  consisted 
of  a  contraction  from  a  12  inch,  18  inch  long  upstream 
reservoir  to  a  2  inch  standard  brass  pipe.  The  reservoir 
for  the  1/2  inch  brass  pipe  was  a  3  foot  section  of  a  3 
inch  brass  pipe.  Both  the  1/2  and  2  inch  brass  pipes 
protruded  1  inch  into  their  respective  reservoirs. 

Figures  (3.2)  and  (3.3)  are  schematic  diagrams  of  the 
upstream  reservoirs  and  indicate  positions  of  the  pres¬ 
sure  taps.  In  addition,  pressure  taps  were  located  along 
the  test  sections  and  their  positions  are  tabulated  in 
Appendix  A,  Table  A-l.  In  order  to  minimize  disturbances 
formed  near  or  on  the  pressure  taps,  the  tap  holes  on 
the  test  sections  were  1/32  inch  in  diameter,  with  each 
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FIGURE  (3. 2) 

UPSTREAM  RESERVIOR  -  2"  BRASS  PIPE 
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FIGURE  (3o3) 

UPSTREAM  RESERVIOR  -  1/2"  BRASS  PI  PE 
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pressure  tap  connected  to  a  purge  line  in  order  to 
remove  trapped  air  from  the  tap  and  from  the  pressure 
lines.  The  pipe  diameter  calibration  in  Table  A-2 
of  Appendix  A ,  resulted  in  the  inner  diameters  being 
0.1721  ft.  and  0.0521  ft.  for  the  2"  and  1/2"  brass 
pipes  respectively. 

3.1.4  PRESSURE  MEASUREMENTS 

Pressure  drops  were  measured  by  two  60  inch 
manometers.  One  manometer  measured  from  0-60  inches 
of  0-dichlorobenzene,  while  the  other  measured  from 
0-60  inches  of  water  resulted  in  a  more  flexible 
measuring  system.  The  temperature-density  relationship 
for  the  0-dichlorobenzene  is  given  in  Table  A- 3  of 
Appendix  A. 

3.1.5  EXPERIMENTAL  FLUIDS 

Two  Newtonian  and  two  polymer  solutions 
different  concentration  were  used  in  this  study.  Water 
and  two  concentrated  aqueous  sugar  solutions  were  used 
to  test  the  experimental  technique.  The  viscoelastic 
solutions  of  0.20%  and  0.01%  having  viscosities  corres¬ 
ponding  roughly  to  that  of  the  Newtonian  solutions  were 
used  as  the  elastic  test  solutions.  This  viscoelastic 
additive,  Separan  AP-30,  is  manufactured  by  the  Dow 
Chemical  Company.  When  the  partially  hydrolized 
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polyacrylamide  (molecular  weight  of  3  x  10  )  is 
dissolved  in  water  the  resulting  solution  usually 
has  a  time  stability:  that  is,  resistant  to  signi¬ 
ficant  degradation  after  an  initial  amount  of  degra¬ 
dation.  This  degradation  is  believed  to  stem  from  a 
partial  or  complete  breakdown  of  the  inter  or  intra 
molecular  forces  or  a  combination  of  these.  However, 

the  0.01%  solution  indicated  a  complete  degradation 

5  -1 

for  shear  rates  exceeding  10  sec  .  This  degradation 
or  decrease  in  drag  reduction  was  found  to  be  irrever¬ 
sible  and  similar  in  magnitude  to  effects  noted  by 
Savins  (65)  with  soap  solutions  that  were  drag  reducing. 

3.1.6  RHEOLOGICAL  PROPERTIES 

Assumed  elastic  and  viscous  properties  of  the 
experimental  fluids  were  estimated  from  data  available 
in  the  literature  (57,66)  and  have  been  reproduced  in 
Appendix  C . 
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3 . 2  PROCEDURE 

3.2.1  PREPARATION  OF  POLYMER  SOLUTION 

Owing  to  degradation  of  the  solutions  it  is 
important  that  each  solution  being  tested  is  treated 
in  a  similar  manner  and  the  following  procedure  for 
dissolving  the  solute  was  adhered  to.  The  solute  was 
added  to  the  tank  filled  to  a  height  of  3  feet  of  tap 
water,  and  stirred  with  an  agitator  for  a  period  of 
several  hours.  Coupled  with  this,  the  pump  was  running 
at  a  moderate  speed  and  continued  to  do  so  for  six 
hours.  Afterwards,  the  solution  was  left  unsheared 
for  the  remainder  of  the  day  and  night. 

The  above  procedure  did  degrade  the  solutions 
to  a  small  extent  however,  significant,  partially 
irreversible,  drag  reduction  loss  for  the  0.01%  solu¬ 
tion  occurred  at  high  shear  rates  in  the  1/2"  pipe. 

This  phenomena  occurred  even  after  only  minutes  of 
taking  data.  Others  (66)  have  noted  that  the  same 
solution  remained  stable  to  further  shearing  for 
several  days. 

3.2.2  DATA  COLLECTION 

/ 

In  order  to  minimize  degradation  pressure  read¬ 
ings  were  taken  in  the  order  of  the  2  and  1/2  inch  pipes 
respectively  and  after  every  few  runs  retakes  were  done 
to  check  if  significant  degradation  had  occurred.  In 
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the  case  of  the  0.2%  solution  the  reruns  checked  within 
experimental  error  except  for  the  0.01%  solution  where 
irreversible,  non-controllable  degradation  occurred 
over  a  short,  continuous  experimental  period  at  high 
shear  rates.  To  investigate  this  phenomena,  experi¬ 
mental  runs  were  made  from  the  low  to  high  shear  rates 
and  the  reverse  employing  a  fresh  solution  at  each  start. 
The  results  yielded  two  different  friction  f actor-Reynolds 
number  curves  over  the  same  range  of  Reynolds  numbers. 
These  effects  are  discussed  in  more  detail  in  Chapter  IV. 
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CHAPTER  IV 

DISCUSSION  OF  RESULTS 

4.1  NEWTONIAN  ANALYSIS 

To  test  the  validity  of  the  experimental  technique 
Newtonian  fluids  were  used  in  both  laminar  and  turbulent 
flow.  As  illustrated  in  Figure  (4.1)  the  friction  factor- 
Reynolds  number  correlations  for  fully-developed  laminar  and 
turbulent  flow  in  a  smooth  pipe  correspond,  within  experi¬ 
mental  error,  to  the  theoretical  values.  The  maximum  abso¬ 
lute  deviation  for  the  experimental  runs  were  within  5%. 

The  results  for  laminar  flow  given  in  detail  in  Table 
(A-2)  of  Appendix  A,  were  used  for  the  determination  of 
the  pipe  diameters.  In  the  entrance  region,  only  four 
pressure  profiles  for  laminar  flow  were  obtained  due  to 
the  inability  to  measure  the  very  small  pressure  drops. 

These  as  well  as  the  turbulent  pressure  profiles  were 
used  to  confirm  the  empirical  laminar  and  turbulent 
contraction  losses  and  entry  length  correlations. 

The  laminar  flow  equation  (1.2.4)  as  shown  on 
Figure  (4.2)  yields  a  straight  line  with  a  slope  of 
32.0  and  an  intercept  equal  to  1.22.  This  intercept 
compares  favorably  with  the  values  of  1.16  and  1.23 
obtained  by  Bogue  (13)  and  La  Nieve  (56)  respectively, 
for  flow  in  capillary  tubes. 
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FRICTION  FACTOR  -  REYNOLDS- 
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As  noted  and  developed  in  the  previous  chapters 
the  relationship  between  the  dimensionless  contraction 
loss  can  in  general,  be  written  as  below  for  Newtonian 
fluids.  Dimensional  analysis  yields  no  indications  as 
to  whether  or  not  the  flow  field  is  laminar  or  turbulent. 
As  a  result,  the  form  of  the  function  A  will  presumably 
be  different  for  the  two  flow  fields. 


— k -  =  A  (  -  )  (4.1.1) 

pV  NRe 

Astarita  (6)  and  La  Nieve(56)  suggest  that  for  high 

Reynolds  numbers  in  tubes  with  or  without  the  presence  of 

a  vena  contracta,  the  dimensionless  contraction  loss  is 

independent  of  the  Reynolds  number  and  remains  constant. 

In  particular,  Astarita  found  that  the  critical  Reynolds 

number  was  146  with  the  proportionality  constant  depending 

upon  the  diameter  ratio  of  the  contraction,  6.  Assuming 

that  this  argument  holds  for  turbulent  flow,  a  plot  of 

log  (AP  )  versus  log  (V)  as  indicated  on  Figure  (4.3) 
c 

resulted  in  a  linear  relationship  having  a  slope  of  2.0 
but  possessing  a  variable  intercept. 

In  Figure  (4.3)  there  is  little  discrepancy  between 
data  from  the  1/2"  and  2"  tubes,  however,  substantial 
differences  in  the  contraction  loss  occur  as  different 


fluids  are  considered. 
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FIGURE  (4„3) 

NEWTONIAN  CONTRACTION 
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2 

□ 

water 
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If  the  density  differences  are  taken  into  account 
the  discrepancies  are  even  larger.  These  facts  are  illus¬ 
trated  on  Table  (4.1)  where  the  intercept  or  Hagenbach  factor 
4>  is  given  for  each  of  the  fluids  in  the  2"  pipe.  Included 
also  in  the  table  are  contraction  losses  and  Reynolds  number 
for  a  velocity  of  10  ft. /sec. 

TABLE  4.1 


VARIATION  OF 

HAGENBACH 

FACTOR 

Solution 

U  (cp) 

Sp . Gr . 

APc (psi) 

*/2 

NRe 

Water 

1.0 

1.0 

1.25 

1.068 

1.6  x  105 

Sugar  #2 

7.0 

1.15 

1.00 

0.920 

26200 

Sugar  #1 

46.5 

1.27 

0.80 

0.830 

4360 

Assuming 

that  the 

correction 

factor  <j> 

is  a  function 

of  the  Reynolds  number,  it  appears  from  the  data  to  be 
directly  proportional  to  the  Reynolds  number.  This  result 
is  similar  to  that  formed  by  Drexler  (28).  The  entrance 
loss  coefficient  as  defined  by  Olson  and  Sparrow  (59) 
yield  values  of  1.136,  0.84  and  0.66  for  water  and  sugar 
solution  #2  and  #1  respectively.  These  values  are  slightly 
larger  than  0.55  obtained  by  Olson  for  a  sharp-edged  contrac¬ 
tion  for  Reynolds  numbers  ranging  from  16000  to  70000. 

The  measured  entry  length  in  turbulent  flow  for 
any  experimental  run  was  less  than  17  diameters.  Latzko's 
theoretical  analysis  reproduced  below  was  developed  for 
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induced  turbulence. 


Lg 

—  =  0.693  N  1/4  (4.1.2) 

D 

The  results  of  this  study  are  in  reasonable 
agreement  with  this  formula  as  indicated  by  Table  (4.2). 

TABLE  4.2 

TURBULENT  ENTRY  LENGTH 

Reynolds  Number  Theoretical  (L^ )  Experiment  (L  ) 

5230  8.5 

37500  9.6 

4.7  x  105  10.20 

2.1  x  105  14.80 

4.2  ANALYSIS  FOR  VISCOELASTIC  FLUIDS 
4.2.1  DRAG  REDUCTION  AND  UNUSUAL  EFFECTS 

Before  proceeding  to  a  discussion  of  entry  lengths 
and  contraction  loss  measurements  for  the  drag  reducing 
material,  a  discussion  concerning  the  method  used  in 
interpreting  the  data  as  well  as  the  unusual  behavior 
under  certain  flow  conditions  is  discussed. 

In  determining  the  friction  factor  as  well  as 
the  shear  stress  at  the  wall  for  each  experimental  run, 
it  was  assumed  that  the  pressure  loss  measured  over  the 
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last  two  pressure  taps  represents  the  developed  pressure 
loss.  Based  on  this,  the  friction  f actor-Reynolds  number 
data  for  each  of  the  solutions  are  presented  on  Figures 
(4.4)  and  (4.5).  The  relationships  as  well  as  indicating 
considerable  drag  reduction  present  unusual  effects  that 
in  the  past  have  not  been  noticed  or  have  not  been  clearly 
defined.  One  of  the  interesting  properties  depicted  on 
Figure  (4.4)  is  that  the  data  from  the  larger  pipe  separate 
from  the  Newtonian  curve  at  a  lower  Reynolds  number. 

Savins  (65)  and  Seyer  (66)  have  noted  this  phenomena  but 
no  explanation  is  available  at  this  time.  Also,  it  is 
clear  that  beyond  a  certain  Reynolds  number  the  friction 
factors  for  the  concentrated  solution  in  the  2"  pipe  as 
shown  on  Figure  (4.4)  have  a  large  amount  of  scatter 
relative  to  the  other  data.  A  more  detailed  analysis 
concerning  this  phenomena  will  be  presented  later  in  this 
section. 

As  the  Reynolds  number  increases,  it  is  noticed 
that  the  data  tend  to  become  parallel  to  the  curve  for 
Newtonian  fluids.  This  may  indicate  that  the  flow  is 
approaching  fully-developed  turbulent  flow  (68) -in  the 
sense  of  an  intermittancy  factor  approaching  unity.  The 
intermittancy  represents  the  percentage  of  time  that  the 
flow  is  turbulent  for  some  given  position  downstream  of 
the  entrance.  The  only  deviation  from  this  approach  to 
a  parallel  relationship  is  best  exhibited  by  Figure  (4.5) 
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FIGURE  (4.4) 

FRICTION  FACTOR  -  REYNOLDS  NUMBER  0.20%  SOLUTION 
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where  the  0.01%  solution  undergoes  a  significant  amount 
of  degradation  beyond  a  Reynolds  number  of  2  x  10  .  This 
degradation  is  reflected  by  a  sharp  decrease  in  the  drag 
reduction  and  the  decrease  continues  until  the  solution 
behaves  essentially  as  a  Newtonian  fluid.  Although  the 
degradation  was  found  to  be  irreversible  over  a  period 
corresponding  to  the  experimental  runs,  it  was  found  to 
be  slightly  recoverable  after  the  solution  had  been 
stored  for  a  prolonged  period.  Savins  (65)  noticed  this 
effect  with  soap  solutions  but  no  mention  has  been  made 
concerning  drag  reducing  systems.  To  comprehend  this 
unusual  effect,  fresh  solutions  were  made  and  data  was 
taken  in  the  order  of  low  to  high  Reynolds  numbers  and 
high  to  low  Reynolds  numbers.  The  results  gave  rise  to 
two  distinct  friction  factor  curves  whose  effects  are 
discussed  relative  to  the  entry  lengths  and  contraction 
losses  later  in  this  chapter. 

To  further  illustrate  the  effects  of  drag  reduction 
and  significant  degradation  reference  is  made  to  Figures 
(4.6)  and  (4.7),  As  indicated  by  the  data  the  0.20% 
solution  for  both  the  1/2"  and  2"  pipes  was  the  better 
drag-reducing  solution.  In  addition,  it  is  clearly 
seen  that  for  the  two  fresh  0.01%  solutions  used  in  the 
1/2"  pipe  the  drag  reduction  of  55%  and  13%  for  the  non- 
degradated  and  degradated  solutions  result  in  a  significant 
distinction  in  the  pressure  losses.  The  highly  degraded 
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solution  resulted  in  only  a  slight  decrease  in  the  pressure 
loss  while  the  non-degraded  solution  resulted  in  the  pres 
sure  loss  being  approximately  equal  to  the  percent  drag 
reduction  as  measured  from  the  friction  factors.  At  low 
velocities  as  indicated  on  Figures  (4.6)  and  (4.7)  drag 
reduction  is  reversed:  that  is,  the  0.01%  solution  is  the 
better  drag  reducing  solution,  due  in  part  to  viscosity 
considerations.  This  is  best  displayed  by  Figure  (4.6) 
where  at  a  velocity  of  5  ft. /sec.  the  pressure  loss 
for  the  0.01%  solution  is  an  order  of  magnitude  less 
than  the  0.20%  solution. 

Typical  data  of  pressure  measurements  along  the 
pipe  are  presented  in  Figure  (4.8)  for  several  flow  rates 
of  the  0.20%  solution  in  the  2"  pipe.  The  unaccounted 
amount  of  scatter  in  the  friction  factor  Reynolds  number 
data  presented  in  Figure  (4.4)  occurred  whenever  the 
flow  rate  exceeded  120  Imp .Gal ./Min.  This  flow  rate 
corresponds  to  the  Reynolds  number  approximately  equal 

4 

to  2  x  10  .  As  indicated  on  Figure  (4.8)  it  is  at 
this  same  flow  rate  where  a  noticeable  pressure  rise 
is  obtained  near  the  entry  of  the  pipe.  The  maximum 
occurs  at  g-  =  8.0  for  all  flows  and  increases  in 
magnitude  for  increasing  flow  rates.  Astarita  (6) 
and  others  (28,59)  have  noticed  this  effect  for  both 
purely  viscous  and  viscoelastic  solutions  and  referred 
to  it  as  a  "vena  contracta  effect" .  For  the  case  of 
laminar  flow  of  viscoelastic  solutions  Metzner,  White 
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and  Denn  (51)  as  discussed  in  Chapter  II,  show  that  the 
stresses  in  rapid  deformation,  short-time  processes  are 
dependent  upon  the  total  deformation.  In  this  limit, 
the  fluid  behaves  like  an  elastic  solid.  For  Newtonian 
fluids  having  small  relaxation  times,  the  Deborah  number 
will  be  small  and  the  solid-like  region  in  the  vicinity 
of  the  contraction  of  cross-sectional  areas  will  not  be 
present.  Metzner  and  White  (47 )  show  at  least  for  laminar 
flow,  that  the  entry  region  can  thus  be  conceived  of  having 
three  distinct  regions:  that  is,  a  region  of  solid-like 
behavior,  a  region  of  boundary  layer  development  and  a 
fully-developed  region.  Under  extreme  conditions,  the 
solid-like  region  could  extend  across  the  entire  duct 
thus  requiring  fluid  fracture  or  slip  of  the  fluid  at  the 
wall.  The  effect  of  this  phenomena  on  the  contraction 
loss  and  entry  lengths  relative  to  the  results  obtained 
in  this  study  is  explained  in  the  following  manner. 

A  qualitative  interpretation  is  made  with 
reference  to  Figure  (4.9).  The  dotted  line  represents 
the  pressure  profile  that  would  exist  if  a  vena  contracta 
was  not  present.  However,  the  viscoelastic  solution  due 
to  its  elastic  behavior,  may  create  a  solid-like  region 
near  the  entry.  This  solid-like  behavior  in  effect  is 
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F I G  U  RE  (4. 9) 

VENA  CONTRACTA  EFFECT 
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equivalent  to  a  reduction  of  the  total  cross-sectional  area 
which  according  to  Bernoulli's  equation  causes  the  fluids 
velocity  to  increase.  This  in  turn  causes  the  pressure 
to  reduce  as  shown.  Once  the  fluid  is  beyond  the  reduced 
cross-sectional  area  it  begins  to  deccelerate  thus  partial¬ 
ly  regaining  its  initial  pressure.  In  this  area  however, 
the  process  is  not  completely  reversible  and  an  extra  pres¬ 
sure  loss  occurs  as  indicated  by  aPLoss*  Beyond  this 
region  the  generally  lower  level  of  acceleration  experienced 
by  the  fluid  allows  normal  development  of  the  boundary  layer 
since  the  response  of  the  fluid  would  be  essentially  that  of 
a  purely  viscous  material.  Thus,  the  overall  effect  is  to 
cause  an  increase  in  the  contraction  loss  by  the  amount 


aLOSS* 

In  this  study,  the  only  "vena  contracta  effect" 
noted  was  in  the  2"  pipe  with  the  0.20%  solution.  No 
effects  were  measured  in  the  1/2"  pipe  with  the  first 
pressure  tap  located  at  L/D  =  5.60. 


4.2.2  CONTRACTION  LOSSES 
(i)  MEASUREMENT 

The  question  arises  as  to  whether  or  not  the 
recorded  pressure  at  the  reservoir  just  above  the  entrance 
is  in  fact  the  actual  pressure  at  the  pipe  entrance. 

Uebler  (77)  and  Astarita  (6)  indicate  that  although  there 

small  pressure  change  due  to  the  height  and  recirculation 
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the  effect  is  negligible  relative  to  the  overall  pressure 
contraction  loss.  Also,  due  to  the  unusual  pressure  curves 
obtained  it  is  necessary  to  indicate  how  the  entry  lengths 
and  contraction  losses  were  determined.  In  addition  to  the 
solid-like  behavior,  pressure  curves  as  shown  on  Figure 
(4.10)  indicate  that  fully-developed  flow  was  not  achieved 
anywhere  in  the  test  section. 

Figure  (4.11)  indicates  qualitatively  how  the 
entry  length  and  contraction  losses  were  measured  for 
both  cases.  Although  curve  1  indicates  the  presence  of 
a  vena  contracta  at  some  point  downstream  the  pressure 
loss  becomes  equal  to  that  of  a  fully-developed  flow 
field.  Extrapolating  this  constant  pressure  loss  per 
linear  foot  back  to  zero  distance  downstream,  a  con¬ 
traction  loss  can  be  measured  as  (AP  ) .  At  the  point 

c 

downstream  where  the  extrapolated  line  separates  from 
the  pressure  profile,  the  distance  measured  is  referred 
to  as  the  entry  length  (L  ).  For  curve  2,  fully-developed 
flow  was  not  achieved  within  the  test  section.  In  this 
case,  the  contraction  loss  was  determined  by  drawing  a 
straight  line  connecting  the  two  lowest  pressure  points 
and  extrapolating  to  zero  distance  downstream  where 
( AP ) 2  is  measured.  To  determine  the  entry  length,  the 
assumed  developed  pressure  loss  over  the  later  portion 
of  the  pipe's  pressure  profile  was  extrapolated  to  beyond 
a  point  where  it  separated  from  the  pressure  profile. 
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The  distance  downstream  where  this  separated  from  the 
pressure  profile  was  taken  as  the  entry  length  (L  )  for 

0  Z 

case  2. 

(ii)  CORRELATIONS 

As  developed  in  Chapter  I  and  employed  in 

Newtonian  analysis  of  turbulence,  the  relationship 

between  the  contraction  loss  and  the  velocity  for 

viscoelastic  solutions  can  be  plotted  in  a  similar 

fashion.  Figures  (4.12)  and  (4.13)  plot  log  (AP  ) 

c 

versus  log  (V)  for  each  concentration  in  the  1/2" 
and  2"  pipes  respectively.  The  results  are  similar 
to  the  previous  analyses  for  turbulent  Newtonian 
fluids  in  that  each  relationship  is  linear.  From 
these  figures  information  related  to  the  vena- 
contracta,  degree  of  elasticity,  and  the  magnitude 
of  the  contractions  relative  to  turbulent  Newtonian 
data  can  be  analyzed. 

One  noticeable  feature  is  that  the  0.01% 
solution  having  a  normal  stress  difference  less  than 
the  0.20%  solution  has  smaller  contraction  losses. 

The  values  of  the  contraction  losses  are  in  fact  approx¬ 
imately  equal  to  the  turbulent  Newtonian  contraction 
loss  for  both  the  2"  and  1/2"  pipes. 


; 


- 


APc(psi.) 


FIGURE  (4. 12) 

CONTRACTION  LOSSES  FOR  1/2"  PIPE 
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FIGUREI4.13) 

CONTRACTION  LOSSES  FOR  2"  PIPE 
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Despite  the  fact  that  the  contraction  losses 

for  the  0.20%  solution  are  larger  than  Newtonian  values, 

the  Figures  indicate  that  they  are  less  than  the  0.01% 

solution  for  laminar  flow  in  capillary  tubes  obtained 

by  Drexler.  Although  not  clearly  indicated  for  the  2" 

pipe  data,  it  is  obvious  from  the  data  for  the  1/2" 

pipe.  Figure  (4.12)  that  the  contraction  losses  approach 

Newtonian  values  for  increasing  velocity.  In  this  case, 

the  relaxation  times  are  in  the  order  of  6.0  x  10  to 
-4 

1.6  x  10  sec.  decreasing  as  flow  rate  increases. 

The  relaxation  times,  measured  in  the  fully-developed 
region,  may  be  indeed  different  from  values  in  the 
entrance  region.  However,  the  data  suggests  that  under 
similar  flow  conditions  the  smaller  the  relaxation  time 
the  contraction  losses  approach  Newtonian  values. 

A  typical  example  is  for  the  1/2"  pipe.  At  a  velocity 

of  60  ft. /sec.  the  0.20%  solution  has  a  relaxation  time 

-4 

equal  to  1.6  x  10  while  the  0.01%  solutions  value  is 
-5 

3.0  x  10  resulting  in  the  contraction  losses  being 
30%  larger  and  equal  to  the  Newtonian  values  respectively. 
One  cannot  extend  this  result  for  a  particular  solution. 
Examining  the  data  for  the  0.01%  solution  in  the  2"  pipe, 
one  finds  the  relaxation  times  having  values  ranging  from 
1.0  x  10  to  3.0  x  10  sec.  yet  have  contraction  losses 
similar  to  Newtonian  values.  It  would  thus  appear  that 
in  general,  the  additive  effects  of  smaller  relaxation 
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times  and  Bernoulli  effect  which  increases  as  the  flow 
rates  increase  cause  the  contraction  losses  for  visco¬ 
elastic  fluids  to  approach  or  become  equal  to  Newtonian 
values . 

For  some  flows,  unusual  vena  contracta  and 
significant  degradation  effects  occurred.  It  was 
previously  indicated  that  the  presence  of  a  vena  contracta 
would  increase  the  contraction  loss.  Although  only  at  a 
flow  rate  beyond  approximately  120  Imp. Gal ./Min.  or  13.5 
ft. /sec.  was  there  any  recorded  effects.  Figure  (4.13) 
indicates  that  there  is  no  sudden  change  in  the  contraction 
loss-velocity  correlation.  If  the  vena  contracta  effects 
commenced  at  this  flow  rate  then  a  noticeable  shift  in  the 
contraction  loss  correlation  should  prevail.  As  this  was 
not  the  case,  it  was  assumed  that  it  existed  closer  to  the 
entry  and  progressively  moved  downstream  as  the  flow  rate 
increased.  Extra  pressure  taps  closer  to  the  entry  are 
required  to  confirm  this  conclusion.  The  significant 
degradation  for  the  0.01%  solution  in  the  1/2"  pipe 
resulted  in  no  measured  difference  to  the  contraction 
loss  as  indicated  by  Figure  (4.12). 

Since  the  solutions  studied  had  weak  to  strong 
viscoelastic  characteristics  it  seems  probable  that  other 
viscoelastic  solutions  such  as  C.M.C.  will  have  contraction 
losses  lying  within  the  results  of  this  study  under  similar 
flow  and  geometrical  conditions.  Considering  the  0.20% 


.rtoi^BXs^ioo  3oX 


ins  srii  oi  :£9aoXo  -qai  o'xusas'iq  s'xixa 


' 


c  nsaz  i.'  .,55  :tei  >r*  -  o  >i  do.-.xv 


■ 

-- 


62 


solution  at  a  velocity  of  10  ft. /sec.  the  data  indicates 
that  there  exists  a  23%  and  74%  increase  in  the  contraction 
losses  for  the  2"  and  1/2"  pipes  respectively.  The  elastic 
effects  thus  create  a  diameter  dependence  not  previously 
noted  for  flow  of  Newtonian  fluids.  This  diameter 
dependence  is  caused  by  the  effect  of  elasticity  being 
more  pronounced  for  smaller  diameter  pipes.  Also,  this 
has  been  measured  (66)  in  terms  of  drag  reduction  by 
noting  that  the  friction  factors  decrease  for  decreasing 
diameters . 

4.2.3  ENTRY  LENGTHS 

At  the  beginning  of  this  chapter  in  Section 
4.2.2  it  was  suggested  that  fully-developed  flow  was 
not  achieved  for  certain  flow  rates  for  the  0.01%  solution 
in  the  2"  pipe.  Despite  this,  the  contraction  losses  and 
entry  lengths  as  discussed  in  more  detail  in  section  (4.2.2) 
can  be  estimated.  The  measured  entry  lengths  are  shown  on 
Figure  (4.14)  and  (4.15)  and  range  from  20  to  180  diameters 
as  the  Reynolds  number  is  increased.  For  this  particular 
solution,  drag  reduction  commenced  at  a  Reynolds  number 
approximately  equal  to  10 ^  and  continued  up  to  a  maximum 
Reynolds  number  of  2.9  x  105. 

As  previously  mentioned,  entry  lengths  in  the  same 
order  of  magnitude  as  those  obtained  by  Rotta  (63)  for 
Newtonian  fluids,  can  in  general  be  predicted  for  visco¬ 
elastic  solutions.  Seyer  (67)  experimentally  predicts 
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long  entry  lengths  to  extend  up  to  Reynolds  numbers  in 
excess  of  30000.  In  this  range,  it  was  also  notic¬ 
ed  by  the  use  of  velocity  profiles  that  the  flow  is 
not  fully-developed.  Also,  Rotta's  data  suggest  that 
beyond  the  transition  region,  that  is,  approaching  fully- 
developed  turbulent  flow  the  entry  lengths  decrease  with 
increasing  Reynolds  numbers.  In  the  laminar  region  the 
entry  lengths  are  directly  proportional  to  the  Reynolds 
number.  For  the  0.01%  solution,  an  increasing  linear 
relationship  as  depicted  by  Figure  (4.14)  results  in 
the  maximum  entry  length  of  200  diameters.  Beyond  this 
no  relationship  can  be  foreseen  due  to  the  limited  amount 
of  data.  It  is  suggestive  however,  that  these  maximum 
entry  lengths  represent  the  entry  lengths  for  a  transi¬ 
tion  region  beyond  which  they  will  decrease  with  Reynolds 
numbers . 

In  addition,  the  two  fresh  0.01%  solutions  used 
in  the  1/2"  pipe  resulted  in  considerable  scatter.  This 
scatter  as  shown  by  Figure  (4.14)  is  noted  for  the  non- 
degradated  solution.  This  non-degraded  solution  had 
a  friction  factor  curve  that  due  to  the  stability  of 
solution  resulted  in  a  separation  from  the  laminar 
Newtonian  curve  rather  than  from  the  turbulent  curve. 
Indicated  by  Figure  (4.14)  this  non-degraded  solution 
has  a  considerable  amount  of  scatter  for  the  low  Reynolds 
number  range  and  appears  to  pass  through  a  maximum  of 
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entry  length  of  60  diameters  occurring  at  a  Reynolds 

4 

number  of  8  x  10  .  Beyond  this,  the  entry  length 
decreases  indicating  an  approach  to  fully-developed 
turbulent  flow.  For  the  degradated  solution,  the  entry 
lengths  continue  to  increase  with  Reynolds  number  until 
a  maximum  entry  length  of  200  is  obtained  at  the  maximum 

5 

Reynolds  number  of  2.6  x  10  .  This  suggests  as  for  the 

2"  pipe  data  that  the  flow  is  in  a  transition  state. 

The  highly  elastic  solution:  that  is,  the  0.20% 

was  the  only  solution  at  least  for  flow  in  the  1/2"  pipe, 

that  resulted  in  a  well-defined  relationship  between  the 

entry  length  and  the  Reynolds  number.  Figure  (4.15) 

indicates  a  definite  maximum  entry  length  of  140  diameters 

4 

occurring  at  a  Reynolds  number  of  2.5  x  10  .  Beyond  this 
Reynolds  number,  the  entry  lengths  decrease  as  the  flow 
field  approaches  fully-developed  turbulent  flow. 

The  presence  of  a  vena  contracta  was  discussed 
in  detail  in  the  preceding  section  and  at  that  time  no 
firm  conclusion  could  be  drawn  relative  to  the  effects 
on  entry  lenths.  It  becomes  apparent  as  shown  by 
Figure  (4.15)  that  the  presence  of  the  vena  contracta 
results  in  considerable  scatter  in  the  entry  lenths. 

This  scatter  results  in  the  entry  lenths  ranging  from 
35  to  130  diameters  for  the  indicated  Reynolds  number 
range.  Despite  this,  the  entry  lengths  are  seen  to  increase 
rapidly  with  Reynolds  number  and  the  maximum  attained 
value  corresponds  to  those  in  the  1/2"  pipe. 
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In  summarizing  all  of  the  above  results,  the 

following  conclusions  are  reached.  First,  the  entry 

lengths  are  always  greater  than  the  corresponding 

Newtonian  values  and  extend  up  to  a  maximum  of  200 

,  5 

diameters  at  a  maximum  Reynolds  number  of  10  .  Second, 
the  results  indicate  that  the  entry  lengths  pass  through 
a  maximum  value  presumably  to  occur  within  the  transition 
region  then,  approach  turbulent  Newtonian  values.  Third, 
the  presence  of  a  vena  contracta  caused  by  the  elastic 
nature  of  the  fluid  and  the  difficulties  in  reproducing 
data  for  the  weak  viscoelastic  solution  in  the  transition 
region  resulted  in  a  significant  amount  of  scatter  in 
the  estimated  entry  lengths. 

4.2.4  DIMENSIONAL  ANALYSIS 


The  concept  of  dimensional  analysis  employed  by 
Astarita  (6)  for  Newtonian  fluids  and  extended  to 
viscoelastic  fluids  (  7  )  in  laminar  flow  is  used  for 
this  analysis.  Considering  the  dimensional  analysis 
as  presented  in  Chapter  II  and  reproduced  here,  one  has: 
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Neglecting  the  effects  of  n':  that  is,  assuming 
it  to  be  a  constant  value  for  each  concentration,  the 
above  equation  can  be  rewritten  as: 
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Astarita's  analysis  considers  only  the  generalized 
Reynolds  number  to  determine  the  Hagenbach  or  Couette 
correction  factors.  Equation  (4.2.1)  indicates  that 
an  additional  parameter,  namely  a  Deborah  number,  may 
be  important.  For  several  values  of  the  Reynolds 
number,  data  showing  the  relationship  between  con¬ 
traction  loss  and  the  Deborah  number  are  presented 
in  Figure  4.16.  For  values  of  the  Deborah  number  less 
than  about  0.1  the  dimensionless  contraction  loss  is 
in  close  agreement  with  values  that  have  been  observed 
for  Newtonian  fluids.  As  the  Deborah  number  approaches 
unity  however,  the  contraction  loss  is  seen  to  increase 
significantly  relative  to  the  low  Deborah  number  or 
Newtonian  value.  At  a  Deborah  number  of  0.60  the 
dimensionless  contraction  loss  for  the  0.2%  solution 
in  the  1/2  inch  pipe  is  1.30  while  the  measured 
Newtonian  value  is  0.83  in  the  same  pipe  at  roughly 
the  same  Reynolds  number.  Equivalently,  these  data 
at  the  noted  conditions,  indicate  a  56%  increase  in 
the  contraction  loss  owing  to  elastic  effects. 

In  Figure  4.17  the  dimensionless  contraction 
loss  has  been  plotted  as  a  function  of  the  generalized 
Reynolds  number.  Considering  the  Newtonian  data, 
plotted  as  the  solid  triangles,  a  distinct  effect  of 
pipe  diameter  is  noted.  In  view  of  the  large  contraction 
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ratio  for  the  system,  which  was  the  same  for  both  pipes, 
the  difference  in  dimensionless  contraction  loss  is 
somewhat  surprising.  Evidently  some  effect  owing  to  dif¬ 
ferences  in  the  entrance  geometry  of  the  two  test  sections 
has  not  been  acounted  for.  In  comparing  data  for  the 
viscoelastic  fluids  with  data  for  Newtonian  fluids  the 
geometry  can  be  accounted  for  by  comparing  data  obtained 
in  the  same  test  section. 

Data  for  the  dilute  solutions  (0.01%)  are  in 
general  agreement  with  the  Newtonian  data.  The  rather 
large  scatter  in  these  data  however,  serve  to  mask  any 
effect  of  geometry  which  may  be  present. 

Comparing  the  0.2%  solution  data,  in  the  1/2  inch 
pipe,  with  the  data  for  the  Newtonian  fluid,  it  is  observed 
that  the  contraction  losses  are  significantly  higher.  This 
effect  is  of  course  eauivalent  to  that  noted  in  Figure  4.16 
at  Deborah  numbers  approaching  unity.  As  a  function  of 
Reynolds  number  contraction  losses  for  the  0.2%  solution 
show  a  significant  decrease  as  the  Reynolds  number  increases. 
This  effect  is  consistent  with  the  observation  that  the 
Deborah  number  decreases  modestly  with  increasing  Reynolds 
number.  For  example  over  the  range  of  Reynolds  numbers 
observed  in  the  1/2  inch  tube  the  Deborah  number  decreases 
from  roughly  0.8  to  a  value  of  0.4  at  a  Reynolds  number  of 


•  ■'  •  *  •  '  i  -  '  - 

- 

.xarfpxri  ^Ijtn*ox5inpxE  sifi  eaaaoi  nox^D.  xinoo  odd  If, rid 


6  9  A 


Comparing  the  data  for  the  0.2%  solution  from  the 
1/2  inch  and  2  inch  pipes,  it  is  clear  that  the  data  from 
the  2  inch  pipe  lie  generally  below  those  from  the  1/2  inch 
pipe.  Equivalently,  the  larger  values  of  the  Deborah 
number  prevalent  in  the  1/2  inch  tube  at  a  given  Reynolds 
number,  serve  to  markedly  increase  the  contraction  loss. 

This  increase  in  contraction  loss  serves  to  more 

than  offset  the  difference  in  contraction  loss  owing  to 

geometry  that  was  observed  for  the  Newtonian  fluid.  Quan- 

4 

titatively,  at  a  Reynolds  number  of  2.2  x  10  the  Deborah 
numbers  in  the  1/2  and  2  inch  pipes  are  0.63  and  0.40 
respectively . 

In  summary,  the  effect  of  elasticity  is  seen  to 
be  a  substantial  increase  in  the  dimensionless  contraction 
loss.  The  available  data,  being  complicated  by  an 
unaccounted  for  geometric  effect,  are  insufficient  to 
quantitatively  define  the  effect  of  elasticity. 
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FIGURE(4.16) 

DIMENSIONLESS  CONTRACTION  LOSS  VS  DEBORAH  NUMBER 
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DIMENSIONLESS  CONTRACTION  LOSS  VS  REYNOLDS  NUMBER 
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CHAPTER  V 

CONCLUSIONS  AND  RECOMMENDATIONS 

5.1  CONCLUSIONS 

For  the  flow  conditions  studied  the  turbulent 
Newtonian  contraction  losses  were  lower  than  those  in 
laminar  flow  and  decreased  in  proportion  to  the  increas¬ 
ing  fluid  density.  This  density  relationship  apparent 
in  the  turbulent  results  do  not  appear  for  the  laminar 
flows.  Entry  lengths  for  the  turbulent  flow  results 
correspond  to  Latzko 1 s  theoretical  formula  with  devia¬ 
tions  occurring  at  high  Reynolds  numbers.  Under  these 
flow  conditions  the  turbulent  entry  lengths  were  less 
than  those  predicted  by  Latzko. 

For  the  0.20%  viscoelastic  solution  in  the  2"  pipe 
a  vena  contracta  effect  was  detected  for  flow  rates  exceeding 
120  Imp.  Gal. /Min.  Although  not  detected  below  this  value 
by  the  pressure  profiles  its  existence  cannot  be  ruled  out. 
The  vena  contraction  is  assumed  to  stem  in  part  from  the 
elastic  nature  of  the  fluid  under  rapid  deformations. 

Under  these  conditions,  the  fluid  resists  the  deformation 
and  as  a  result  sets  up  a  solid-like  response  near  the 
entry  of  the  tube.  The  overall  effect  is  an  increase  in 
the  contraction  losses  relative  to  turbulent  Newtonian 
data.  As  the  flow  rate  increases,  the  discrepancy 
between  the  viscoelastic  and  the  Newtonian  data 
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diminish.  This  is  due  to  the  increased  Bernoulli 
effect  which  in  the  case  for  the  0.01%  results  in 
contraction  losses  equal  to  the  Newtonian  values. 

The  entry  lengths,  although  suffering  from 

significant  scatter  show  a  definite  increase  from 

those  expected  for  Newtonian  fluids.  Entry  lengths 

in  the  order  of  200  diameters  were  found  for  the  0.20% 

and  0.01%  solutions  at  Reynolds  numbers  approximately 

4  5 

equal  to  2  x  10  and  3  x  10  respectively.  The  maximum 
attained  value  presumably  occurs  within  the  transition 
region.  Beyond  the  maximum,  the  entry  lengths  at  least 
for  the  0.20%  solution,  approach  turbulent  Newtonian 
values . 

The  use  of  dimensional  analysis  indicate  that 
the  dimensionless  contraction  losses  depends  upon  a 
dimensionless  group  which  includes  the  relaxation  time 
of  the  material,  6.  It  has  been  found  however,  that 
this  single  dimensionless  group  as  insufficient  to 
quantitatively  interpret  the  results.  In  addition 
to  the  above  group,  the  Reynolds  number,  although 
depicting  the  results  for  laminar  flow,  was  not 
capable  of  describing  the  turbulent  results .  It  is 
apparent  that  other  dimensionless  groups  such  as  the 
elastio-inertial  number  must  be  considered  in  more 
detail  before  a  definitive  analysis  of  the  entry 
phenomena  can  be  obtained. 
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5 . 2  RECOMMENDATIONS 

The  following  recommendations  are: 

1.  A  modification  should  be  made  to  the  pressure  taps 

on  the  pipes  near  the  entry  to  enable  one  to  confirm  that 
the  measured  pressure  in  the  reservoir  above  the 
entry  can  be  considered  equal  to  that  in  the  pipe. 

Also,  with  more  taps  located  near  the  entry,  a  more 
detailed  analysis  can  be  obtained  concerning  the 
"vena  contracta  effect"  or  the  solid-like  behavior 
that  was  experimentally  present  for  the  0.20% 
solution  in  the  2"  pipe.  These  taps  would  have  to 
be  less  than  an  ^  =  8.0.  As  an  added  feature, 
the  pressure  measurements  would  indicate  the 
position  where  the  boundary  layer  changes  from 
a  laminar  to  a  turbulent  condition  as  detected 
by  an  inflection  point  on  the  pressure  curves. 

2.  Experimental  studies  should  be  conducted  with  other 
moderately  or  highly  elastic  solutions  so  as  to 
yield  a  more  comprehensive  analysis. 

3.  Longer  test  sections  especially  for  the  larger 
diameter  pipes  is  required  in  order  to  determine 
one  way  or  another  if  the  assumption  that  the  pres¬ 
sure  drop  over  the  last  section  does  indeed  represent 
the  developed  pressure  loss. 
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4.  Velocity  profiles  would  aid  in  the  determination 
of  the  boundary  layer  growth  and  the  evaluation 
of  the  shear  stress  at  the  wall  in  the  entry 
region.  These  would  allow  one  to  theoretically 
determine  the  magnitude  of  the  normal  stresses 
and  the  previously  mentioned  boundary  layer.  In 
addition  a  more  exact  picture  would  be  available 
concerning  the  solid-like  behavior  in  the  entry 
region  for  high  flow  rates  in  small  or  large  pipes. 

5,  Coupled  with  the  above,  the  relationship  between  the 
normal  stress  differences  with  radial  position 
would  enable  one  to  employ  other  semi=empirieel 
approaches . 
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EQUIPMENT  DIMENSIONS 

With  reference  to  Figure  (A-l) ,  Table  (A-l) 
specifies  the  downstream  and  upstream  pressure  tap  loca¬ 
tions  for  both  the  2"  and  1/2"  smooth  brass  pipes.  The 
total  length  for  the  2"  and  1/2"  test  sections  of  42.5' 
and  16.5'  respectively  resulted  in  the  last  pressure  tap 
being  37  and  45  diameters  upstream  from  the  end. 


2" 

BRASS 

1/2" 

BRASS 

REFERENCE 

TAP  LOCATION* 
(FEET) 

L/D 

TAP  LOCATION* 
(FEET) 

L/D 

1A 

-1.000 

5.81 

-1.000 

19.20 

2A 

-0.500 

2.91 

-0.500 

9.60 

3A 

0.000 

0.00 

-0.250 

4.80 

4A 

- 

- 

0.000 

0.00 

1 

0.500 

2.91 

0.292 

5.61 

2 

1.535 

8 . 93 

0.451 

8.66 

3 

3.084 

17.93 

0.920 

17.70 

4 

6.084 

35.40 

1.829 

35.20 

5 

16.084 

93.50 

4.860 

93.50 

6 

26.084 

151.80 

7.927 

152.50 

7 

36.084 

210.00 

10.919 

210.00 

8 

_ 

_ 

14.095 

271.00 

* 


Reference:  Entrance  to  brass  pipe  test  section. 
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TABLE  2 
TEST  SECTION 

Using  sugar  solution  number  1  in  laminar  flow 
the  calibration  of  the  pipe  diameters  was  done  in  the 
following  manner.  Using  the  manufacturer's  data  for  the 
2"  and  1/2"  brass  pipes  and  knowing  the  flow  rate  an  expe¬ 
rimental  friction  factor  is  determined.  This  experimental 
value  when  compared  to  the  theoretical  value  resulted  in 
the  absolute  deviation  of  5.0  and  2.9%  for  the  1/2"  and  2" 


pipes  respectively. 

Thus ,  the 

manufacturer ' s 

values  of 

0.1721'  and 

0.0521' 

for  the  I. 

D.  of  the  pipes 

were  used 

in  this  study.  Table  (A—  2 )  lists  the  experimental  and 

theoretical 

friction 

factors . 

PIPE  (INCHES) 

NRe 

f 

exp . 

■^theoretical 

DEVIATION 

1/2 

326 

.  0491 

.  0494 

+.0003 

1/2 

470 

.  0341 

.0348 

+.0007 

1/2 

538 

.  0297 

.  0282 

- . 0015 

1/2 

860 

.0184 

.  0180 

+.0005 

1/2 

1500 

.  0107 

.  0095 

-.0012 

1/2 

1900 

.  0084 

.  0088 

+.0004 

2 

416 

.  0385 

.0396 

+.0011 

2 

607 

.  0264 

.  0257 

-.0007 

2 

882 

.  0182 

.0184 

+.0002 

2 

1164 

.  0138 

.0142 

+.0004 

2 

1302 

.0123 

.  0128 

+.0005 

2 

1995 

.  0080 

.  0081 

+.0001 

j 
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TABLE  A- 3:  MANOMETER  FLUID  DATA 

Densities  for  the  manometer  fluids  were  obtained 
from  Seyer  (66) .  The  fluids  used  were  mercury  and  O-dichloro- 
benzene.  The  temperature  density  relationships  for  water 
and  O-dichlorobenzene  are  given  on  Figure  (A-2) .  The  density 

3 

of  mercury  was  assumed  to  be  a  constant  value  of  13.550  lbm/ft  . 
TABLE  A- 4:  SUGAR  SOLUTION  DATA 

For  each  of  the  two  sugar  solutions,  a  temperature 
viscosity  and  a  temperature  specific  gravity  relationship 
was  determined  by  use  of  a  Brookfield  Viscometer  and  a 
Westfall  balance  respectively.  Since  most  of  the  experi¬ 
mental  runs  were  between  a  room  temperature  of  65.0°F  to 
85.0°F  the  relationships  were  determined  within  the  same 
range  as  shown  on  Figures  (A-3)  and  (A—  4 )  for  the  sugar 
solution  1  and  2  respectively. 
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Temperature  (  F) 


60. 


64. 


76. 


80 
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1.278 


1.277 


1.276 


1.275 


1.274 


1.273 


1.272 


1.271 


1.270 


Temperature  (°F) 


Specific  Gravity 
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7.4 


7.2 


7.0 


P-t 


6.8 


6 , 6 


6.4 


6.2 


6.0 


Temperature  (  F) 


60 


65.0 


80,0 


Viscosity  (C. 
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CALIBRATION  DATA 
>  J  - 

FLOW  METER  CALIBRATION 

T  - -  1  ~ 

Calibration  from  0  -  65.0  Imp.  Gal. /Min.  was 
accomplished  by  direct  weighing  and  from  30.0  to  200.0 
Imp.  Gal/Min.  by  one  employing  the  Newtonian  Von  Karman 
correlation 


j  =  4.0  log  (NRe/F )  -  .40 

Calibration  from  weighing  and  by  the  correlation 
were  within  1.0%  and  1^%  respectively  of  the  full  scale 
repeatability  of  the  flowmeter. 

For  flow  rates  ranging  from  0-25  Imp.  Gal. /Min. 
the  converter  was  set  at  span  number  2.  This  high 
sensitivity  setting  resulted  in  a  more  accurate  reading 
for  these  low  flow  rates.  The  results  tabulated  on 
Table  (B-l)  are  plotted  on  Figure  (B-l) .  For  the  higher 
flow  rates  the  converter  setting  was  altered  to  span 
number  1.  Both  weighing  and  employing  the  Von  Karman 
equation  as  tabulated  on  Tables  (B-2)  and  (B-3)  are 
plotted  on  Figure  (B-2) .  In  this  case,  the  weighing 
served  to  check  the  accuracy  of  the  pressure  drop 
readincrs  . 
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TABLE  B-l 


0-25  IMP 

.  G.P.M.  CALIBRATION 

(HIGH 

SENSITITY) 

Settings : 

Amplifier 

-  IV/ cm. 

Room  Temp.  72 

Converter 

-  Span  No 

.  2 

Flow  Temp.  70 

RUN  NO. 

WEIGHT 

(LB) 

TIME  FLOW 

(SEC)  (IMP.G.P. 

RECORDER 

M.)  DIVISIONS 

1 

40.0 

58.70 

4.08 

3.30 

2 

60 . 0 

54.65 

6 .58 

5.30 

3 

60 . 0 

48.20 

7.47 

5 . 80 

4 

75.0 

58.00 

7.75 

6 . 00 

5 

75 . 0 

35.10 

12 . 81 

10.05 

6 

100.0 

45.80 

13.10 

10.50 

7 

100.0 

43.20 

13 . 89 

11.00 

8 

90 . 0 

40.00 

13.50 

10.75 

9 

100.0 

27.00 

24.45 

19.25 

■A 
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TABLE  B-2 

0765  IMP.  G.P.M.  CALIBRATION  (LOW  SENSITIVITY) 

Settings:  Amplifier  -  1  volt/ cm.  Room  Temp.  72°F 

Converter  -  Span  No.  1.  Flow  Temp.  70°F 


RUN  NO. 

WEIGHT 

(LB) 

TIME 

(SEC) 

1 

55.00 

44 .00 

2 

75.00 

37 . 30 

3 

60.00 

29 . 60 

4 

65.00 

32.60 

5 

90.00 

44.30 

6 

80.00 

25 . 00 

7 

140.00 

25.10 

8 

140 . 00 

19 . 00 

9 

150.00 

21.10 

10 

150.00 

25.80 

11 

190.00 

17.20 

FLOW  RECORDER 

(IMP. G.P.M.)  DIVISIONS 


5.50 

2.00 

12.08 

3.25 

12.15 

3.25 

11.98 

3.20 

12.20 

3.30 

19.20 

5.00 

33.50 

8.70 

44.25 

11.35 

42.60 

11.20 

58.10 

15.00 

66.30 

17.20 

O.vX 
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TABLE  B-3 


0-170  IMP.  G.P.M.  CALIBRATION 

Settings : 

Amplifier  -  IV/cm. 

Room  Temp: 

7  2°F 

Converter  -  Span  No.  1 

Flow  Temp: 

7  0°F 

RUN  NO. 

PRESSURE  DROP/20  FEET 
(IN.  OF  WATER) 

FLOW 

(IMP.  G.P.M.) 

RECORDER 

DIVISIONS 

1 

110.90 

165.30 

44.80 

2 

91.50 

146.20 

40.00 

3 

80.50 

136.80 

37.50 

4 

67.70 

124.00 

34.00 

5 

61.10 

117.00 

30.40 

6 

45 . 80 

99.80 

26.30 

7 

32.70 

82 . 80 

22.20 

8 

23.30 

68.50 

18.40 

9 

17.60 

58.65 

15.80 

10 

13.00 

49.50 

13.10 

11 

9 .40 

41.25 

11.00 

12 

5.60 

30.90 

8.00 

13 

3.53 

23.80 

6.00 

14 

6.13 

32.45 

8.00 

15 

10.25 

43.30 

11.00 

16 

17.50 

58.30 

14.75 

IMP.  GAL. /MIN. 
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FIGURE  B-l 

Low  Flow  Rate  Calibration 

Amplifier:  IV/cm.  Room  Temp:  72°F 

Converter:  Span  No.  2  Flow  Temp:  70°F 
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High  Flow  Rate  vs.  Recorder  Divisions 
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PHYSICAL  PROPERTIES 

TABLE  C-l;  PHYSICAL  AND  RHEOLOGICAL  PROPERTIES 

The  flow  curves  used  for  the  two  concentrations 
of  Separan  AP-30  were  obtained  from  Seyer  (66)  and  are 
reproduced  in  Figure  (C-l)  .  The  densities  used  were  the 
densities  of  water. 

The  relaxation  time  as  a  function  of  the  shear 
stress  at  the  wall  for  the  0.20%  and  0.01%  Separan  AP-30 
were  obtained  from  the  available  literature  of  Oliver  (57) 
and  Seyer  (66)  and  as  plotted  on  Figure  (C-2) .  As  noted, 
for  the  0.01%  solution,  there  is  a  discrepancy  between 
Oliver's  and  Seyer 's  data  which  is  accounted  for  by 
degradation.  Since  the  0.01%  solution  is  this  study 
resulted  in  two  discrete  friction  factor  relationships, 
it  is  to  be  expected  that  each  would  possess  a  different 
relaxation  time  -  shear  stress  curve;  however,  no  capil¬ 
lary  jet  thrust  apparatus  was  available  and  as  a  result 
Seyer 1 s  relationship  was  chosen  to  represent  both  the 
degradated  and  non-degradated  solutions. 


,  -.3  .tulos  #10.0  srfd’  10I 


FIGURE  (C-l) 


C-2 


in 

o 

rH 


o 


ro 

o 

i — l 


<N 

O 


M. 

Ii 


( • jsd) 


8V/D  (Sec. 


C-3 


(  *  O0S)  9 


Figure  C-2:  Relaxation  Time-Shear  Stress  Data 


Appendix  D 


Newtonian  Pressure  Profiles 


■  M 


D-l 


Newtonian  Pressure  Profiles 

The  tabulated  experimental  pressure  drops 
given  on  Tables  (D-2)  to  (D-6)  are  plotted  on  Figures 
(D-l)  to  (D-8).  From  these  pressure  profiles  one  is 
able  to  measure  entry  lengths  and  contraction  losses 
and  calculate  the  Reynolds  numbers  and  friction  factors. 
These  quantities  are  tabulated  on  Table  (D-l) . 

The  following  outlines  the  order  in  which  the 
figures  appear. 
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VISCOELASTIC  PRESSURE  PROFILES 

Due  to  significant  amounts  of  degradation  and 
the  unusual  types  of  pressure  profiles  for  the  solutions 
studied  new  solutions  were  used  extensively.  This  pro¬ 
vided  a  check  for  reproducibility  and  clearly  indicated 
the  above  effects  on  entry  length  and  contraction  loss. 

Table  (E-l)  lists  the  order  of  the  figures  and 
describes  when  fresh  solutions  were  used.  Table  (E-6) 
to  (E-9)  tabulates  the  experimental  pressure  drops  that 
are  plotted  in  Figures  (E-l)  to  (E-18).  From  these 
figures  one  is  able  to  measure  entry  lengths  and  contrac¬ 
tion  losses  as  well  as  calculating  Reynolds  numbers  and 
friction  factors.  These  measured  and  calculated  variables 
are  tabulated  on  Tables  (E-2)  to  (E-5)  inclusive. 


{  3\7iin  X9  V>  -1  If  JX  J  X  -  1 
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E-2 


FIGURE  NO 

E-l  to  E- 

E-4  to  E- 

E-8  to  E- 

E-ll  to  E 

E-14  to  E 


TABLE  (E-l) 


ET-597  PIPE  DIA. 

%  CONC.  (IN) 


COMMENTS 


3  0.20 

7  0.20 

10  0.01 


-13  0.01 


-16  0.01 


2  -  Fresh  solution 

-  Low  to  high  flow  rates 

-  Figure  E-3  reproduce 
data  with  another  fresh 
solution  (Runs  10  to  20). 

1/2  -  Used  solution  from  above 

-  Low  to  high  flow  rates 

2  -  Fresh  solution 

-  Low  to  high  flow  rates 
for  runs  (1-14) . 

-  Runs  (15-22)  reproduce  the 
above  data,  with  another 
fresh  solution.  Located 
in  Chapter  4,  Figure  (4.10) 

1/2  -  Fresh  solution 

-  Low  to  high  flow  rates 
for  runs  1  to  11 
inclusive . 

1/2  -  Fresh  solution 

-  High  to  low  flow  rates 
for  runs  12  to  24 
inclusive . 

1/2  -  Fresh  non-degradated 

solution 

-  Low  flow  rates  only  for 
runs  25  to  29  inclusive. 


E-17  to  E-18 


0.01 


:  O.  *  X: 
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E-3 


TABLE  (E- 2) 


%  Cone.:  0.20  Fluid  Temp.:  68.0°F 

Pipe  Dia.  :  2  in.  Room  Temp.:  70.0°F 


RUN  NO. 

ve  ,  .2 

—  X  it) 

Re 

fxlO3 

APC 

(IN.H20) 

Le/D 

1 

35.93 

1236 

18.033 

2.32 

8.70 

2 

39.09 

2659 

8.874 

5.55 

17.50 

3 

50 . 89 

5326 

4.014 

16.05 

93.00 

4 

54.30 

9564 

2.316 

36.50 

40.60 

5 

51.56 

14364 

1.710 

54.00 

75.70 

6 

42.12 

20117 

1.513 

82.50 

35.00 

7 

49.80 

25225 

1.085 

Jr 

116.00 

122.00 

8 

47.59 

19211 

1.422 

85.00 

110.00 

9 

53.21 

29643 

0  .900 

154.00 

134.00 

10 

46.99 

19490 

1.422 

89.50 

93.00 

11 

38.47 

24786 

1.365 

119.50 

139.00 

12 

29.47 

18306 

2.183 

78.00 

35.00 

13 

30.58 

31797 

1.337 

149.00 

35.00 

14 

33.04 

25658 

1.503 

110 . 00 

47.00 

15 

30.78 

19826 

1.983 

76.50 

35.00 

16 

36.46 

2267 

1.072 

5.30 

47.00 

17 

38.91 

1265 

1.661 

2.30 

8.70 

18 

51.17 

17302 

1.466 

73.50 

64.00 

19 

39.14 

23297 

1.423 

106.00 

116.00 

20 

38.37 

31068 

1.133 

152.50 

116.00 

■ 

, 

dv 

E-4 


TABLE 

(E-3) 

%  Cone . : 

0.20 

Fluid  Temp . : 

(Xi 

<sO 

• 

o 

O 

•n 

Pipe  Dia.: 

1/2  in. 

Room  Temp . : 

En 

0 

o 

• 

o 

RUN  NO. 

ve  ,  n2 

—  x  io 

N  ' 

Re 

fxlO3 

APC 

(in.h2o) 

L  /D 
e 

1 

89 . 06 

3948 

4.706 

34.00 

40.30 

2 

71.79 

7872 

3.091 

80.00 

86.50 

3 

70.13 

14411 

1.873 

168.00 

93.00 

4 

62.92 

21641 

1.441 

293.00 

92.00 

5 

6 

53.56 

32053 

1.170 

508.00 

110.00 

7 

51.16 

39108 

1.020 

673.00 

124.00 

8 

n 

43.01 

45799 

1.029 

813.00 

110.00 

y 

10 

38.56 

41266 

1.238 

24.75 

86.30 

11 

34.59 

50910 

1.115 

37.00 

62.50 

12 

30.65 

57913 

1.123 

38.75 

67.20 

13 

29.60 

71683 

0.962 

46.50 

24.00 

14 

25.50 

82836 

0.968 

66.75 

28.80 

15 

24.31 

88178 

0.952 

74.40 

25.00 

16 

21.30 

92145 

1.027 

/ 

73.25 

13.50 

17 

24.30 

81525 

1.019 

61.75 

30.60 

18 

26.47 

73679 

1.435 

47.50 

44.00 
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TABLE  (E-4 ) 

%  Cone.:  0.01  Fluid  Temp.:  65.5°F 

Pipe  Dia.:  2  in.  Room  Temp.:  73.0°F 


RUN  NO. 

ve  ,  „2 

—  X  10 

N  . 

Re 

f  xlO  3 

APC 

(in.h2o) 

L  /D 
e 

1 

1.30 

33841 

5.145 

2.30 

35.00 

2 

1.38 

45755 

4.433 

3.60 

40.60 

3 

1.44 

70790 

3.902 

7.70 

23.20 

4 

1.40 

58013 

4.183 

5.40 

46.50 

5 

1.47 

103941 

3.547 

15.50 

35.00 

6 

1.52 

104632 

3.334 

15.20 

24.20 

7 

1.56 

123797 

3.113 

24.00 

46.50 

8 

1.62 

138128 

2.875 

31.00 

26.10 

9 

1.54 

172660 

3.078 

46.50 

175.00 

10 

1.50 

207192 

3.158 

62.00 

185.00 

11 

1.54 

239134 

2.945 

88.00 

93.00 

12 

1.57 

290587 

2.760 

128.00 

174.00 

13 

1.81 

203394 

2.206 

62.00 

163.00 

14 

1.40 

288343 

3.424 

130.00 

174.00 

15 

1.31 

30215 

5.140 

1.20 

2.90 

16 

1.75 

174387 

2.400 

43.00 

151.00 

17 

1.81 

208919 

2.210 

69.00 

174.00 

18 

1.85 

221523 

2.104 

78.00 

163.00 

19 

1.83 

241724 

2.124 

87.00 

180.00 

20 

1.86 

270213 

2.019 

112.00 

182.00 

21 

1.84 

296112 

2.057 

136.00 

188.00 

22 

1.81 

290587 

2.112 

132.00 

186.00 

. 

, 

IS- .  I 

TABLE  (E-5) 


%  Cone.:  0.01  Fluid  Temp.:  67.5°F 

Pipe  Dia.:  1/2  in.  Room  Temp.:  74.0°F 


RUN  NO. 

N  ' 

Re 

fxlO3 

APC 

(IN.H20) 

Le/D 

1* * 

4.91 

28517 

3.523 

12.00 

48.00 

2* 

5.94 

55893 

2.274 

43.00 

58.00 

3* 

6.00 

77566 

2.142 

82.00 

58.00 

4* 

6.15 

107794 

1.960 

175.00 

48.00 

5* 

5.74 

145437 

2.153 

325.00 

25.00 

6* 

5.67 

169391 

2.160 

430.00 

24.00 

7 

5.43 

193916 

2.304 

57.00 

9.60 

8 

5.18 

221292 

2.481 

76.00 

9.60 

9 

5.02 

250950 

2.582 

95.00 

5.00 

10 

4.79 

274905 

2.793 

119.00 

5.00 

11** 

4.48 

288023 

3.151 

131.00 

19 . 00 

12 

4.98 

292015 

2.569 

1370.00 

153.00 

13** 

4.74 

275475 

2.839 

1220.00 

210.00 

14** 

4.42 

246388 

3.279 

1055.00 

173.00 

15** 

4.39 

232129 

3.352 

927.00 

180.00 

16** 

4.33 

196768 

3.521 

615.00 

170.00 

17** 

4.16 

174524 

3.848 

470.00 

38.00 

18** 

4.51 

157985 

3.345 

328.00 

27.00 

19** 

4.46 

138023 

3.472 

275.00 

9.60 

20** 

4.29 

102091 

3.867 

153.00 

9.50 

21** 

4.18 

71863 

4.260 

75.00 

14.30 

22** 

3.86 

34220 

5.406 

17.50 

5.50 

23** 

3.92 

28517 

5.360 

8.50 

5.50 

24** 

3.76 

51330 

5.411 

43.00 

7.70 

25* 

6.18 

37072 

2.214 

20.00 

» 

19.30 

26* 

5.62 

27376 

2.747 

9.80 

14.50 

27* 

5.83 

13688 

2.796 

1.55 

86.00 

28* 

5.61 

7699 

3.227 

0.75 

96.00 

29* 

4.94 

19962 

3.646 

5.00 

5.50 

*  These  runs  were  used  to  determine  the  friction  factors  for 
the  non-degradated  solution  on  Figure  (4.5). 

*  These  runs  were  used  to  determine  the  friction  factors 

for  the  degradated  solution  on  Figure  (4.5). 
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PRESSURE  DROPS  (in.  WATER) 
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FIGURE(E-2) 

Pressure  Profiles  for 
0.20%  Separan  A  P-30  in  2"  Pipe 
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FIGURE(E-IO) 

Pressure  Profiles  for 

0.01%  Separan  AP-30  in  2"  Pipe 
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FI  G  URE(E- 15) 

Pressure  Profiles  for 

0.01%  Separan  A  P-30  in  1/2"  Pipe 
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Appendix  F 
Sample  Calculations 
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The  sample  calculations  given  below  are  based  on 
Run  1  of  the  0.20%  Separan  AP30  data  in  the  2"  pipe. 


Table  G-l: 

For  Run  1 : 

Q  =  18.00  Imp.  Gal. /Min. 

AP  =  3.34  in.  water. 

L  =  10  ft. 

For  the  2"  pipe  the  bulk  velocity  is  given  by: 

(18.00)  (4.0)(  .  1336  )  (1.2) 

V  = - 5 - 

(60.0)  (tt)  ( .  17  21 )  Z 


=  2.069  ft. /sec. 


The  shear  stress  at  the  wall  is  obtained  by 

converting  the  pressure  drop  to  psf  by 

(3.34)  (14.7)  (144.0) 

AP  =  - 

(33.91)  (12.00) 

=  17.375' (psf. ) 

and 

(17.375)  (0.1721) 

T  = - 

W  (4. 0)U0. 0) 

=  0.0748  (psf.) 

Now  the  friction  factor  can  be  determined  as 

(0.0748)  (2.0)  (3  2. 2) 

f  =  j 

(62.4)  (2.069; 


0.01803 


F-2 


With  the  value  of  the  shear  stress  at  the  wall  and 
use  of  Figure  (C-l)  in  Appendix  C  the  flow  behavior  index, 
consistency  index  and  the  laminar  shear  rate  can  be  determined. 


n'  =  0.557 

K'  =  0.004216 


8V 

D 


96.0 


The  generalized  Reynolds  number  can  now  be  calculated 


as 


Nd  ' 
Re 


(.1721) *557  (2.069)  2,0  *  557  ( 62 . 4 )  81,0  *557 
(0.004216)  (32.2) 


=  1236.51 


From  Figure  (C-2)  in  Appendix  C  the  relaxation  time 
of  0.02990  sec.  is  obtained.  This  is  now  used  to  determine 
the  normal  stress  difference 


(P-,-P0<5)t  =  (0.0299)  (2.0)  (0.0748)  (175.0) 

11  2  2  W 

=  0.7796  (psf.) 

The  Deborah  number  and  inertia-viscous  number 
can  now  be  determined  as 


N 


DEB 


V  9 
D 

(2.069)  (0.0299) 
(0.1721) 


(0.359) 


F-3 


(P. , -P00)  g 
11  22  w^c 

PV2 

(.7796)  (32.2) 

(62.4)  (2.069) 2 

0.094 

To  determine  the  dimensionless  contraction  loss  and 
entry  length  one  must  obtain  these  values  from  Figure  (E-l) 
Appendix  E . 

(AP  )TT„  =  2.32  (in.  water) 
c  V 

L  =  1.50  ft. 

e 


and 


N 


IV 


Thus,  the  dimensionless  contraction  loss  becomes 


(APc}VE 


pV2/ q 


(2.32)  (14.7)  (144.0)  (32.2) 
(12.0)  (33.91)  (62.4)  (2.069)2 


=  1.458 

and 

L  1.50 

e 

D  0.1721 


8.70 


Appendix  G 
Nomenclature 


' 
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a 

constant 

A,  A' 

arbitrary  functions 

Ar  Ai 

arbitrary  functions 

b 

distance  between  parallel  flat  plates 

c(c,f) 

constant  defined  by  (16) 

D 

pipe  diameter 

d.  . 

ID 

deformation  rate  tensor  defined  as 

h  U  .  v.  +  a  .v. ) 
i  D  D  1 

f 

Panning  friction  factor 

F  . 
entry 

frictional  loss  in  the  entry 

Fdev 

frictional  loss  in  developed  flow 

F  ,  frictional  loss  in  the  upstream  reservoir 

upstream  ^ 

F„  frictional  loss  due  to  elasticity 


AE 

elastic  stored  energy  per  unit  mass 

5C 

r  ' 

j  i  i  *~7  a  f  t  •  lbm 

constant  32.174  - - - 

sec . ^ lbf  j 

K 

constant  defined  by  equation  (1.2.1) 

K' 

consistency  index  defined  by  equation  (2.2.1) 

K2 

Hagenbach  correction  factor  defined  by  equation 
(1.2=6) 

Kl 

%e 

Couette  correction  factor 

AKE 

kinetic  energy  change 

L 

distance  downstream  from  the  contraction 

Le 

entry  length 

Le 
( — ) 
VD;ve 

viscoelastic  dimensionless  entry  length 

' 


" 


G-2 


Le 

(^>pV 

purely  viscous  dimensionless  entry  length 

n' 

flow  behavior  index 

Deb 

dimensionless  Deborah  number  defined  as  (— ) 

D 

N 

ws 

dimensionless  Weissenberg  number  defined  as 
(pll"p22)w 

Tw 

NRe 

dimensionless  Reynolds  number 

NRe 

dimensionless  generalized  Reynolds  number 

P 

hydrostatic  pressure 

P.  . 

ID 

deviatoric  stress  component  defined  by  equation 
(2.2.3) 

AP 

pressure  loss  in  developed  flow 

iPDEV 

pressure  loss  in  developed  flow 

aploss 

defined  by  Figure  (4.9) 

AP 

TOT 

total  pressure  drop  measured  from  the  contraction 

AP 

c 

contraction  loss 

R 

pipe  radius 

s 

normal  stress  power  law  index  defined  as 

(Pll_P22^w  =  const-  (-$-) 

S  .  . 

D 1 

strain  tensor  (39) 

T 

w 

shear  stress  at  the  wall 

T.  . 

ID 

component  of  stress  tensor 

u 

instantaneous  axial  velocity 

+ 

u 

dimensionless  velocity  defined  as  (u/u*) 

u* 

fxictional  velocity  defined  as  \/Tw/p 

G-3 


V 

average  velocity 

y 

distance  from  pipe  wall 

+ 

y 

dimensionless  distance  defined  as  (yu*/v) 

Greek  Symbols 


6 

dimensionless  diameter  ratio 

aloss 

defined  by  Figure  (4.9) 

5  .  . 

ID 

P 

Kronecker  delta 

contravariant  convective  derivative  (52) 

density 

0 

relaxation  time  defined  by  equation  (2.4.7) 

4 

Hagenbach  correction  factor 

V 

kinematic  viscosity 

a(n,s)  constant  defined  by  equation  (1.2.2) 

y  viscosity 


A 
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